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A major portion of the present-day energy requiremenxs 
is met by hydroelectric and fossil fired thermal power stations 
However, the increasing demand for power and the current oil 
crisis make it imperative to explore other sources of energy* 
Nuclear fission energy seems to he a promising alternative* 
Efficient operation of existing nuclear power plants and dev- 
elopment of plants for future needs requires an understanding 
of the nuclear reactor and other auxiliary systems.. In this 
work, modelling of a pressurised heavy water reactor li- s 0 cl J* 'X' 
out. The model is then simulated on a digital computer to 
study its dynamic - uncontrolled and controlled - hehaviour* 

A reactor of the CANDU type - CANDU is a heaxyr water 
moderated* pressurised heavy water cooled natural uroii&4.m. 
fuelled press'ure tube reactor - is studied in this thesis. Ttie 
system for producing sjfeeam comprises U-type shell and tube heat 
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exchangers? where primary coolant’ transfers heat to water and 
a steam drum. One leg of the shell of the heat exchanger 
forms the riser for natural circulation unit while the other 
forms part of the feed water circuit. A mathematical model 
represented hy a set of coupled nonlinear, ordinary, differ- 
ential equations is formulated for the above mentioned system - 
reactor, heat exchanger and steam drum with natural circulation 
loop - from physical principles. The system is then simulated 
on a digital computer to study its dynamic behaviour when 
operating at full load, it is acted upon by disturbances like 
step change in the area of the steam valve, step change in the 
temperature of the feed water, etc. In order to design 
controllers so that the average temperature of the coolant 
in the reactor, the pressure of the primary coolant circuit, 
and the water level in the steam dium are maintained at 
prescribed values, the nonlinear model is linearised by a 
first order Taylor series expansion around full load conditions. 
The resulting linear, coupled, ordinary differential equations 
with constant coefficients are then expressed in state space 
form to facilitate the application of optimal control theory. 
Uncontrolled responses of the nonlinear model are compared 
with those of the linear model for equivalent disturbances to 
validate the linearization procedure. 



!Ih9 problem of selection of controllers is then 


formulated as a linear regulator problem, where an objective 
function involving a quadratic function of states and control 
inputs is minimised subject to system’s governing equations. 
Firstly a state feedback controller of proportional type - 
control inputs are linear, constant combinations of all states 
of the system ~ which drives the system states to prescribed 
values depending on magnitude of disturbance, is obtained. 
However, this is not suitable if the disturbances^ are not 
Tuaeasurable or if the parameters of the system vary from the 
nominal values. To overcome these drawbacks, a proportional 
and integral type of controller which forces selected states 
to zero in the presence of constant disturbances is obtained'. 
Such a controller requires a knowledge of all state variables. 
Since some of the state variables are not available for control 
purposes, a simpler control scheme based on available states 
is desirable. To achieve «uch a controller, the linear 
regulator problem is treated as a function minimization prob- 
lem in the parameter space of the feedbacks, for which the 
necessary condition of optimality is that the gradient of the 
objective function with respect to each allowed feedback 
must vanish. Using explicit expressions fqr the value of the 
objective function and its gradient with respect to each 
allowable feedback, a descent method based on the gradients 
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is used to minimize the objective function and to achieve 
optimality. This approach leads to the design of controller 
of specified or preselected configuration and it is utilised 
to design analog controllers which minimize a performance 
index. Such a controller is implemented on the linear model. 
Controlled responses of the system with the simpler, easy- 
to-implement control scheme compare favourably with those 
where control . scheme requires feedback of all the states, 
lastly, the optimal analog controller obtained from the 
analysis of linear system is implemented on the nonlinear model 
and the controlled responses of the nonlinear model are 
provided for comparison. 





CHAPTER 1 
INTRODUCTION 

A major p/artion of the present day energy requirements 
is met by hydroelectric and fossil fuelled thermal power stations 
However increasing demand for power and the current oil crisis 
make it imperative to explore other forms of energy such as 
nuclear energy and solar energy. Nuclear fission energy seems 
to he a promising source for future electrical power generation. 
In the U.S.A., U.S.S.R., Canada, U.K, and other industrialised 
nations, a number of nuclear power plants are operating success- 
fully. In India, two nuclear power stations are in operation 
and tw* more are being coiiimissioned. The number of fsuch 
installations is bound to in^crease rapidly in the foreseeable 
future with many countries opting for nuclear power, ^ 

1.1 ROLE OP COMPUTERS IN NUCLEAR POWER PLANTS 

The expected increase in the number of nuclear power 
plants with a consequent increase in their sizes will place 
great emphasis on the development of instrumentation system 
that will increase plant availability and prevent outages. 
Punctionally , the instrumentation provides data which can be 
used for fuel management, detection and correction of p»wer 
distribution changes in the core, control of plant parameters 
like st0eam pressure, temperature, water-level in the steam 
drum etc. In most plants the above functions are accomplished 
by analog devices. However it has been felt that a coinputer 
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could be readily applied to the nuclear plant to perform 
many of the functions itemized above^, 

i'cdoy proness-control computers are being increasingly 

used in the operation of nuclear power plants. Both start-up 

and shutdown routines and steady state operation are being 

entrusted to computers in some power stations . In addition, 

a plant computer can relieve the operator of the tedious job 

of logging data and can simultaneously analyse the data if an 

operational problem exists. Of the 56 US thermal stations 

under constriiction in 1970, 48 had computers that performed 

digital monitoring, 53 bad computers that calculated performance 

data and 14 bad computers for starting tbe turbine-^. Tbe 

success achieved in conventional power stations has led to 

tbe use of computers in nuclear power plants, the trend being 

to utilise the computer to assimilate data, make calculations 

and as an operational aid. Its use in control and safety 

functions is generally avoided'^. However power plants like 

Douglas Point Generating Station have used computers for direct 

digital control. Ihe promising results obtained from Douglas 

Point G.S. have led to the application of computer systems in 

Ontario Hydro’s Pickering G.S. and Bruce G.S. and Gent illy 

5 

Nuclear G.S. in Canada . When operating experience- justifies 
their reliability, economy, and accuracy, nuclear plant computers 
will be used for various automatic control functions. In 
anticipation of this, efforts are being made to apply advanced 



3 


technology to the plant control system and st-udies^ have 
been conducted for optimal digital control of nuclear reactors. 

1.2 iitbbA®^RE review 

Analysis and design of optim'al feedback control systems 

require a detailed understanding of the dynamical behaviour of 

7 

the power station. The work of Chien et al, is directed towards 
this objective. Chien et al. developed a mathematical model for 
a drum type natural circulation boiler. Dynamic behaviour of 
the boiler was described about a steady state operating level 
by perturbing the nonlinear mass, momentum and heat balance 
equations and the heat transfer equations. Uncontrolled 
transient responses were obtained for various inputs. Conven- 
tional analog controllers for closed loop control were selected 
based on a comparison of controlled responses of the plant 
with uncontrolled responses for equivalent distvirbances. The 
work of' Chien et al. initiated the use of digital computer 
for studying the dynamical behaviour of power plants. 

8 

Nicholson , following an analysis similar to that of 
Chien et al. obtained open loop responses of an oil-fired 
boiler. The order of the mathematical model was reduced 
analytically by neglecting the isolated higher order modes. 
Uncontrolled responses of the original model and the reduced 
model were compared. Using the reduced model, optimal digital 
control was obtained based on the minimisation of a performance 
index related to system*s state variables and system inputs. 
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The dynamic behaviour of the system with optimal control was 
compared with the behaviour of the system with conventional 
analog control. In a subsequent work^, turbo-alternator was 
included to study tjie combined boiler- turbo-alternator system. 
Performance of the system with optimal digital control and 
with conventional control was then investigated, 

4nderson^° studied the dynamic behaviour of a natural 

circulation boiler. Transient responses using conventional 

analog control, direct digital control and optimal digital 

control were obtained with a view to study the possibility of 

improving the performance of the plant by using an optimal control 

11 

system. A similar study was performed by Anderson et al. for 
a 200 MW boiler. In this analysis, random input distiurbances 
and measurement noise were taken into account. Thompson and 
Adams et al^^ developed mathematical models for drum type 
boiler and once-through boiler respectively and the responses 
of the model agreed with the responses obtained from field 
tests. An IBM report^^ details the simulation models for 
hydro, thermal and nuclear power systems. 

Dynamic analysis of a nuclear boiler based on a steam 
generating, pressuxe tube, heavy water moderated design was 
performed by M'Pherson et al^^ A detailed mathematical model to 
represent plant dynamics was described with speoinl reference^^^^^^^^^^^^^^^ 
to the boiling chanhel and steam drum dynamics. A simplified 
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model. was derived for analogue computer s'tudies* Effect of 
analog control systems on the performance of the plant was 
then studied. 

Optimal control theory has heen used to study the problems 

associated with the reactor such as optimal xenon shut do'vm^^^ 

17 

xenon spatial oscillations , reactor start up and transient 

18 

response, lipinski provides a detailed survey of the works- 
done in optimal control of the reactor. 

19 

Kalley suggested four applications of dynamic prog~ 
ramming teclinique to nuclear reactor- These are: optimisation 
of poison distribution, optimisation of overall plant efficiency, 
design of optimal control programs and determination of flow 
distribution through a heat exchanger. Kalley outlined the 
optimal solution to a minimum energy start-up problem. 

20 

Ash used dynamic programming to derive a functional 
eq\iation which would cause a boiling water reactor to be driven 
to its equilibrium condition in minimum time by continuously 
moving the control rods. 

21 

Using Pontryagin’ s principle, Mohler obtained the 

optimal open loop control to study the minimum time control of 

22 

neutron density, Sekhar and Weaver investigated optimal 
closed loop control: of a reactor by using a set of equations 
linearised around the nominal trajectory alohgwith the minimisation 
of a quadratic performance index. 
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Dwicorabe defined a control variable which, was equal 
to the product of flux density and reactivity and hence made 
the neutron kinetics equations linear. The performance index 
included a term of reactivity times flux squared. Optimal 
closed loop control was obtained by using the maximum principle 
for nuclear reactor load control, lipinskit studied the 
optimal digital control of nuclear reactors with estimation of 
state variables. In the above studies, only the reactor 
kinetics with suitable approximations such as prompt jump, 
and one- group delayed neutron model was considered. 

PA 

Ciechanowicz*^^^ investigated tvo types of digital control 
schemes for the linear model of the Halden Boiling Water Reactor 
one involves the control strategy calculation for the overall 
dynamic system and in the other, the overall system was spli"^ 
into t^o systems characterised by smaller niimber of variables, 
the interaction between the two sub-sys'^ems being included by 
the use of cross-coupling controllers. The splitting produces 
the advantage of reduced computational effort. Dynamic 
behaviour of HB¥R under different modes of control was compared 
with the behaviour of the reactor without control. 

6 

Bjirlo et al. studied the transient response of the 
pressure in a boiling heavy water reactor. Experimental results 
which compare digital control and existing analog control for 
equivalent perturbations were provided. It was then concluded 
that digital control handles the perturbations effectively. 
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25 

Ciechanowicz’ s work is a combination of the earlier 

C O A 

studies ^ ^ for the investigation of multi-level control of a 
power plant. A 500 MW BWR power plant, treated as three 
coupled cores was studied for controlled and uncontrolled 
behaviour. Overall control strategy was investigated for 
different sampling periods of digital control and for the 
case where constraints were imposed on reactivity controller 
signals. Further, the influence of noise upon the overall 
dynamics of the power reactor was studied. 

26 

SLivinsky and Weaver'^ used the state variable 
technique to design non-interacting controllers for a system 
with m- inputs and m-outputs and illustrated the design of 
controllers for coupled core reactors where each core is 
controlled by the controller corresponding to that particular 
core. 

Beneznai and Sinha'^ ’ investigated the design of 
digital controllers with provision to adopt to changes in 
plant parameters. (Kie response of the reactor was approximated 
by that of a low order model in a desired manner, the model • 
parameters being chosen to minimize the deviation between 
the system and the model responses. Optimal feedback 
parameters for the model were then computed by means of the 
matrix jRiccati equation. 
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1.3 OQTIIKE Aim SCOPE OE THE THESIS 

In India, there are two nuclear power stations; the 
one at Tarapur has a "boiling water rear tor and the other in 
Rajasthan has a pressurised water reactor. A few plants 
similar to the one in Rajasthan are expected to he commissioned 
in the near future. The research described in this thesis was 
undertaken with the objective of applying modern control theory 
to design a control system for a pressurised water reactor. 

Chapter 1 is devoted to a review of the application of 
computers in power plants and of the optimal con+rol of 
nuclear reactors. 

A brief description of a pressurised heavy water reactor 
is made in Chapter 2. In Chapter 3 the nonlinear ordinary 
differential equations which describe the' dynamic behaviour 
of the plant are derived from physical principles. Uncontrolled 
transient responses for selected disturbances are also obtained 
by simulating the nonlinear model on an IBM 7044 digital computer. 

In Chapter 4, the nonlinear equations developed in the 
previous chapter are modified and are linearised around the 
steady state operating point at full power. This leads to the 
representation of the system by a set of linear, ordinary 
differential equations. These equations are manipulated into^^^^^^^ 
state space form to yield a linear time invariant system for 
the reactor. Transient responses are then obtained for 
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di star "banc es identical to those assumed for the nonlinear 
model. A comparison of the uncontrolled responses of the 
two models is made to validate the linearisation procedure. 

In Chapter 5? the control problem is formulated as 
a linear regulator problem with the objective function as a 
quadratic function in the estate and control variables. 
Controllers of the proportional type and of the proportional 
and integral type are obtained. Behaviour of the plant with 
these controllers providing the feedback control are discussed 
when the plant is acted upon by constant disturbances, 

The controller obtained in Chapter 5 assumes that 
all state variables are measurable and hence are available 
for the' construction of the control law. This may not be 
possible in many physical situations. In such cases, a 
controller of preselected configuration is desirable. 

Therefore in Chapter 6, a controller whose structure is 
specified is obtained. The design of this controller is 
formulated as a function minimisation problem in the parameter 
space of the feedbacks. A descent technique based on the 
gradient of the performance index with respect to each allow- 
able feedback is used to design the controller. Transient 
behaviour of the plant with this control is compared with 
that obtained in Chapter 5 for equivalent disturbances. 

Further, this controller is suitably modified and is 
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incorporated in the nonlinear model developed in iOdapter 5 
and the transient responses of the plant are discussed for 
constant disturbances* 

Finally in. XJhapter 7 the results of the present 
work are summarised and suggestions for further work are 
provided* 
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DESCRIPTION OP PRESSURISED HEAVY WATER REACTOR 

2.1 INTRODUCTION 

A conventional thermal power station and a nuclear 
power station differ mainly in the heat generating unit - 
the former utilising a furnace where chemical changes of the 
fuel (coal or oil) produce heat energy and the latter using a 
nuclear reactor where nuclear changes in the fuel liberate the 
heat. The hea.t thus liberated is used to produce steam in 
a steam generator. Most of the other equipments such as steam 
turbines and electric generators are quite alike for both the 
power plants. In this chapter brief descriptions of nuclear 
power reactors and more specifically a pressurised heavy water 
reactor of the CANDU type are presented. 

2.2 DESCRIPTION OP NTJCLE/iR REACTOR 

The basic features of a nuclear reactor are the reactor 

core in which heat is produced, a coolant system to remove the 

heat from the core and a heavy shield surrounding the core to 

contain all emanations that might otherwise be harmful to 

29 

human life or health . The core consists of a fuel which 
contains uranium, plutonium or thorium atoms. When a neutron 
is absorbed by a fuel nucleus, there is a finite probability 
of splitting this nucleus into fragments resulting in the 
release of two or more neutrons. This process is known as 
fission. The neutrons produced due to fission may take 
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part, in addition to causing further fission, in several 
reactions; they may he absorbed in core material other 
than the fuel or may be lost by leaking out of the core. 

For the fission process to be sustained, for each fuel 
nucleus capturing a neutron and undergoing fission, a 
minimum of one neutron, on the average, must be produced 
and absorbed by another fuel nucleus to cause fission. Since 
a large amount of energy (about 2C0 million electron volts or 
0.32x10 Watt sec. per fission) is released in fission, 
a sustained fission process provides a steady heat source. 

A salient feature of the fission process is that all 

neutrons are not relaased simultaneously. Most of the neutren*^ 

“12 

are emitted within about 10 seconds after fission. These 
are known as prompt neutrons. However, a small portion, less 
than one percent, of the total neutrons is given off an 
appreciable time after the fission has occurred and these 
neutrons are called delayed neutrons. These are emitted due 
to the radioactive decay of certain fission fragments, known 
as precursors. The half life of these precursors vary from 
0.05 second to 55-6 seconds. The delayed neutrons are of 
tremendous importance since these facilitate the control 
of chain reaction, which otherwise would have been quite 
difficult, if not impossible. 

Depending on the average speed of the neutrons in 
the reactor, reactors are classified as fast, intermediate 
and thermal reactors. Heutrons released in fission are at 
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high, energies (i.e, have high velocities). Neutron energy- 
on the average is greater than one million electron volt. 
When such neutrons cause most e;f the fissions, then such 
reactors are called as fast reactors. However, neutrons 
having energies less than about 0.1 electron volt have greater 
probability of causing fission than fast neutrons. In some 
reactors the fast neutrons are slowed down to thermal energies 
(0.025 electron volt) by collisions with certain materials, 
known as moderator, before they are absorbed by the fuel. 

Such reactors, where thermal neutrons maintain fission, are 
called thermal reactors and the fuel-moderator arrangement 
is known as core. 

Heat generated in the core is removed by circulation 
of a coolant. Depending on the type of fuel, moderator and 
coolant , reactors are classified into various categories such 
as the boiling water reactor, pressurised water reactor, 
gas cooled reactor, etc. 

A large number of commercial nuclear electrical power 

generating systems utilise water cooled and moderated thermal 
30 

reactors . Iwo basic types are generally used - pressurised 
water reactors which have no b-ulk boiling of coolant in the 
core and boiling type which permits bulk boiling in the core. 
Of the two operatitig atomic power stations in India, the ' 
one at Tarapur uses a boiling water reactor whereas the 
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reactor of Rajasthan Atomic Power Station (RAPS) is of the 
pressurised variety. Systems, similar to RAPS, are under 
construction at Madras and at Rarora. 

2.3 DESCRIPTION OP AN ATOMIC POWER STATION USING A 
PRESSURISED WATER REACTOR 

31 

An outline of the power station is as follows : 

The power station consists of a CANDU I Cana d ian 
Deutrium Uranium) type reactor supplying steam to a turbo- 
generator of 200 MW generating capacity. CANDU is a heavy 
water moderated, pressurised heavy water cooled, natural 
uraniiim fuelled pressure tube reactor. The reactor has a 
horizontal cylindrical tank and tube assembly. This assembly 
is known as calendria. Natural uranium in the form of 
pressed and sintered pellets is shes-thed in Zircaloy tubes 
to form cylindrical fuel elements. Nineteen such elements 
are assembled to form a fuel bundle. The calendria ■ tube 
contains thick walled coolant tube, having a fuel bundle 
and primary coolant, viz., heavy water at high pressure. The 
shell of the calendria contains heavy water at low pressure. 
This heavy water is used as moderator. 

The two functions of the heavy water namely as a 
moderator and as a coolant are separate, each having its 
own closed circulating system- The fuel coolant system, 
called the primary heat transport system, is a high pressure 
and high temperature circuit with no bulk boiling. The 
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moderator circuit is a low temperature, low pressure circuit. 
These two arc insulated from each other in the reactor core 
by a gas annulus between the calendria tube and the coolant 
tube. 

Primary heat transport system transfers the heat 
generated in the fuel elements to heat exchangers where the 
heat is utilised to generate steam in a secondary coolant 
system. The primary system, shown in Pig. 2.1, consists essen-- 
tially of two identical loops in series, one at each end of the 
reactor. Each loop includes four boiler assemblies in parallel 
and four pumps in parallel. Such an arrangement circulates 
the primary coolant, through alternate coolant tube assemblies, 
in opposite directions. Bach boiler assembly consists of ten 
vertical, U type, tubes and shell heat exchangers in parallel 
and a steam drum. 

The hot heavy water from the reactor enters the tube 
sides of the heat exchangers at the bottom of their recir- 
culating sections after a transport delay of 5.6 seconds. 

Then it returns from the bottom of the preheater legs of the 
heat exchanger to the reactor. The transport delay associated 
with the flow from the steam genera,tor to the reactor is 
of the order of 7.2 seconds. 

The secondary coolant system, shown in Pig. 2. 2 consists 
of boiler feed water entering the shell side of the each heat 
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exchanger at the hottOEi of the preheater leg* The feed water 
is heated to or near the saturation temperature in the 
reduced diameter length of the leg. It boils in the upper 
portion of the leg and enters the steam drum above it. Here 
water and steam are separated. The water at saturation 
temperature flows through down-comers to the bottom of the 
heat exchanger's recirculating leg where it boils and returns 
to the steam drum. At fioll load, design pressure at the 
steam drum is 569 psig. 

Saturated steam (0,-24 percent maximum wet) from the 
steam drum is supplied to the turbine unit at 565 psig. The 
turbine has one sin^e flow high pressure element a.nd tw^o 
double flow low pressure elements. The exhaust from the low 
pressure unit is condensed in a surface condenser. The 
condensate after six stages of heating, returns to the preheat 
leg of the heat exchanger, thereby completing the secondary 
coolant circuit. 

The quantity of steam produced in the secondary 
system varies with reactor power. Power is controlled by 
adjusting the reactivity. This can be accomplished by 
changing the moderator level or by adjusting the neutron 
absorbing rods, called control rods, or by adding or removing 
a poison to the moderator. In the present study, however, 
it is assumed that reactivity changes are brought about by 
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movement of control rod only. 

Pressure of the heavy- water coolant at reactor outlet 
is maintained at 1250 psia by feeding or bleeding heavywater 
to or from the primary heat transport system. Water level 
in the steam drum is controlled by adjusting the feedwater 
flow. The level of water is designed to vary linearly with 
load. But in this work, the water level in the steam drum 
is controlled to be at the centre of the drum. 

Discharge valves are provided in the steam mains to 
release steam to dump condenser in case of pressure rise in 
the drum. This facility is provided in case the turbine is 
unable to accept enough steam. In addition to the above 
mentioned controls, the turbine and the reactor have additiona,! 
controls for the safe and smooth operation of the plant. 

In this thesis, a section of the power plant consisting 
of the reactor, the primary heat transport system, the heat 
exchangers, the steam drum and steam valve in the steam mains 
is simxfLated. A mathematical model for the above system is 
formulated, the details of which are given in the next chapter. 



CHAPTER 3 
NOH-LIEEAR MODEL 


3.1 INTRODUCTION 

In the previous chapter, different units of a 
nuclear power station based on a pressurised heavy water 
reactor were briefly described. Since the power station has 
to function as a system, rather than as a collection of 
individual pieces of equipmeniSf it is essaatial that the 
system behaviour as a whole be reasonably well described 
and understood. This is required in order to establish the 
necessary control strategies so that the performance of the 
system meets the specifications imposed from operational, 
safety and other consideration. In other words, the aim is 
to design a system i.e., the process plus control, that x\rill 
do a prescribed job in the ’best’ possible way, the best in 
one case may be determined by economics and in another, by 
transient responses to a set of expected disturbances. In 
this chapter, the dynamic behaviour of the pressurised heavy 
water reactor is studied and a mathematical model is obtained 
as first step towards the design of a control system. 

3.2 NEED EOR A MODEL 

The first step towards achieving this objective is to 

7 

completely understand the process being controlled . This 
begins with a phenomenological description of the process by 
following the nuclear, physical and chemical changes that 
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occur along the process path and ends when the time dependent 
functional relations between the input (controlling) and output 
(controlled) variables are obtained. This results in a math- 
ematical description of the system, the solution of which 
predicts the dynamic behaviour of the system to various distur- 
bances. Unfortunately, an accurate mathematical description 
of the entire system is quite difficult, if not impossible to 
achieve, since the system is complex and contains numerous 
variables which are unwieldly to manipulate. Nonlinearities 
and uncertainities in certain physical ‘phenomena aggravate the 
complexity of the problem. 

To overcome the above difficulty, it is necessary to 
make some simplifying assumptions. These assumptions may be 
classified under two categories depending on the nature of 
simplifications; 

(a) Certain physical phenomena are too complex to admit an 
exact mathematical description. Hence some semi- empirical or 
approximate formulation is required, examples being heat transfer 
phenomenon in boiling and two phase flow. 

(b) The second type of assumptions stem from the need to 
obtain a solution Ttithin a reasonable time for the set of 
equations which describe the system. Equations for conservation 
of momentum and energy and for neutron population in the reactor 
are in general nonlinear partial differential equations, the 
exact solutions being limited to simplified cases only. Hence 
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simplifications may be made in tbe eq^uations tbemselves by 
eliminating nonlinearities or by eliminating partial differ- 
entiation through appropriate techniques such as finite 
difference methods. In the present work, these simplifications 
along with certain additional assumptions are combined such 
that a set of ordinary differential equations describe the 
dynamic behaviour of the plant. 

3.3 DYN/iMIC BEHAVIOUR 

Eor the purposes of developing mathematical expressions 
which describe the system, the power plant is divided into 
the foIlDwing sub-systems; 

i) the reactor and the fuel element with the coolant 
flowing around it, ■ 

ii) the heat exchanger which is further divided into the 

recirculating leg, the boiling leg and the preheat leg, 

iii) the steam drixm with the natural circulation loop, 

iv) the valve at the steam mains, and 

v) the primary heat transport system 

Bellowing sections describe the formulation of the 
dynamic behaviour of the above sub-systems. 

It has generally been assumed that the coolant in the 
primary heat transport system is in the liquid phase onlyi 
Changes in its density with temperature variations are 
neglected, except while calculating the primary circuit 
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pressure. Also for the t¥0~phase (steam-water) flow, 
specific volume of the mixture v is calculated as 
V = XVg + (1-x) v^ 

where 


X 




= quality of the mixture 

= specific volumes of saturated liquid and 
vapour respectively 


The quality x is calculated as the mass fraction of steam in 
the two-phase mixture. 


3«3»1 Reactor 

Dynamic behaviour of, the reactor can be characterised 
by neutron population which is dependent on reactivity. Control 
mechanism is the primary source for providing reactivity varia- 
tions. Changes in the temperatures of the fuel and of the 
coolant in the reactor however introduce reactivity, which is 
known as feed back reactivity. Thus the temperatures of the 
coolant and the fuel affect indirectly the behaviour of the 
reactor. 

Temperature of the fuel depends on the power at which 
the reactor is operating and on the temperature of the coolant 
in the reactor. It can be visualised that the coolant temper- 
ature at inlet to the reactor core is determined by the system 
external to the reactor, i.e., heat exchanger and steam drum. 
Hence, changes in operating conditions and variations in the 
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performance of the external system influence the operation 
of the reactor. 


For purposes of control studies, it is adequate to 


represejnt the reactor by a simple point model with six groups 

,30 


of delayed neutrons' 


Bie point model kinetics equations are 


uiil = n(t) + 2 - c.(t) 


and 

where 


dc^(t) 


i=l 


^i 


n(t) - Cj_(t), i = 1,2, ...6 


n(t) 

6k(t) 

^i 


neutron density 
= reactivity 

: total delayed neutron fraction 
: neutron life time 

: decay constant of the i"^^ neutron precursor 
: concentration of delayed neutron precursor 
of group i 

■ delayed neutron fraction of group i 
If neutron density and concentration of the precursor 
are normalised with respect to their initial steady state 
values (at time t=0), the equations take the form 


Ci(t) 


P. 


n' (t) + E 

4 “ ^ \ / 


dt 


i=l 


( 3 . 1 ) 


dc! (t) 

= Xj_ [n» (t) - c'j^(t)], i = 1,2,. -.6 


( 3 . 2 ) 
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where n’ and cj are normalised values, with initial conditions 
1 

being = 1.0 for i = 1,2, ...6 and n' =1,0. 

Ihe number of neutrons in the core is proportional 
to the n\mber of fissions occuring and for 33dL0^^ fissions per 
second, one watt of thermal power is produced. The power output 
of a reactor is then proportional to the number of neutrons in 
the core at a given instant of time and so the symbol n’ (or n) 
is used to designate the neutron level, with the implication 
that a power level is involved. 


Heat generated in the fuel due to fl'tisions occurring 
there is transferred by conduction in the fuel and the clad. 
Heat transfer from fuel to clad is a. mixture of convection and 
conduction. Coolant flowing around the clad removes the heat 
by convective transfer. 


Conduction of heat in a solid medium has been discussed 
52 55 

in the literature’'^ * , With internal heat generation, the 

temperature distribution in the solid is given by the 
following equation: 


3T ^ 

"ST pc 


pc 


where 

p = density of the medium 
c = specific heat of the medium 
■ k = thermal conductivity of the medium 
q = heat energy developed in unit volume and time 
= Laplace operator 
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Thermal conductivity k is in general a function of temperature. 
The above nonlinear partial differential equation is solved 
by finite difference technique. With the steady state 
temperature distribution thus obtained, the fuel dynamics is 
approximated and the average temperature of the fuel is 
represented as 


^ = 107.6( m ). 0.1938 (3.3) 

^100 

■where n’ ^ refers to the full design value normalised with 
100 

respect to the initial steady value. Appendix A details the 
calculation of the coefficients from steady state conditions. 


Considering the coolant in the channel as a single 
lumped system represented by the average temperature of the 
coolant between inlet and outlet of the core and assuming that 
the heat loss by the coolant tubes to the moderator is 
negligible, the average temperatuire of the coolant can be 
expressed as 

= 5.44 (87/2723 ) + 2-30 (1^1-1^^) 

The coefficients are calculated ba.sed on steady state conditions, 

Rea-ctivity 6k(t) in Bq. (3*l) can be divided into two 

parts i.e., one due to feed back other due to 

control mechanism 6k . 6k-, is written as 

c fb 

6k^^(t) ^ ATg^^(t) + AT^j(t) , (3.5) 
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where 

a = temperature coefficient of reactivity 
for the coolant 

= temperature coefficient of reactivity for 
the fuel 

and 

AT , AT = changes in fuel and coolant temperatures 

£ll Q# O 

from their respective steady state values. 

Heat Exchanger 

Eor the purpose of analysis, the heat exchanger is 
divided into three units - recirculating leg, "boiling leg and 
preheat leg .and each unit is simulated separately to study 
its dynamical "behaviour. It is also a.ssumed that all the 
eighty heat exchangers which are connected to 8 steam drums 
are identical in their performance. 

(a) Recirculating leg 

Firstly, the recirculating leg is approximated by a 
set of n lumped sections. Mass and energy balances for the 
three media viz., the primary coolant, the tube material 
and the secondary coolant (water-steam mixture) yield the 
steady state temperature distribution. Based on the steady 
state distribution, the riser leg is approximated as a single 
Itimped system. Appendix B outlii^es the calctilation of the 
dynamic behaviour of the single lumped system and a comparison 
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of the responses obtained for the single lumped system with , 
that of the distributed system approximated by n lumped 
section. 

The recirculating leg can therefore be represented by 

T^ = 0.45 (V^w2^ 

T^2= 5.296 (12-^2^ ~ 0.15152x10“^ exp (p/225 ) (5}^2“^sat)'^ 
h^ 20. 854x1 0“5 exp (p/225 )(T^2“^sat)'^ + 

Wd^^f'^r^ ^ Pr) (3.6) 

Similar analyses are carried out for the boiling leg 
and the preheat leg (Appendix B) and for these legs the 
dynamic equations can be written as 


(b) Boiling leg 

Tj = 0.577 (Ij-Tj) - 1.743 

= 2.702 - 0.1243x 10'5 exp(p/225)(l„3-Igat)'^ 

A = '{ 12.81x10"^ exp (p/225) 

> ^0-102/ (s-T) 

(c) Preheat leg 



1.76 (Tj-T^) - 0.3204 (l4-'^w4) 

f 

0.6558 (T^-i;^) - 5.95 (T^^- ■ 


^sat+\i 


(203.5 (T 


^sat+^mi- 


“^mi(V^V^'l^^sat“^nii^))>/ * (3.8) 


The symbols have the same meaning as those in Appendix B. 
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3 •5*3 Steam Drum and Natural Circulation Loop 

There are eight steam drums in the plant. It is 
assumed tha.t these eight steam drums behave identically 
i.e,, the flow rates into and out of these drums are 
identical and the liquid is at the same level in a,ll these 
drums. Following additional assumptions are made: 


i) There is no spatial variation of temperature in the 
drum and the temperatures of liquid and vapour phases 
are assumed to be equal to the saturation temperature 
corresponding to the instantaneous pressure, 

ii) The steam separation in the drum is perfect, 

iii) Bj^bble formation in the liquid is neglected, 

iv) Heat capacities of drum internals and drum material 
are neglected. 


Based on these assumptions# equations for the 
conservation of mass and energy are manipulated and the 
dynamic behs-viour of the drum pressure and water level 
can be written as 


P 


where 


A’ 

B» 


C 


= (A‘-«B'-C»-D'b2h^)/(E»+P>b^+G’b^)c (3-9) 

^g ^fg 

“ ^’^fg \ “ ^fg ^f '*■ "^f ^fg^^^mi'*’\ ^b ^4 

= (hy-h^) v^g b^ hj,) 
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i)' 

E’ 

E" 

pt 


V h. „ 7 0 

3 ^ fg r 

dv 

(;?r^)(li- -v^ B») 

g dp fg fg ^ 


M 


/(M 


V ) 


I dp f g dp g 

\ - ^fg Ilf + Hfg) \ % 


2 

r 


t (V^f ) '"fg + i^g} \ p 

h dF 

^fg^V 

4 (-fg/^fg) - 5 ! 


'fg 


M 


g 


= (ird^/S-yd)?- P 


g 


and --=(ttd ^S+yd ) ilp^ 


♦ 

y —(A'+B'-C'r-E* )/(E’xp») 


where 


” ’'mo 

’'mo “ ’'ml ■’■ ^b» t \ + S j) 
= ' = ’'r 

’'r =’'a + (’’a ^ \ p) 

=''a'^ig 

Vg.E» 


(5.10) 
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E” = 
F» = 


M 


dhj 


dv. 


■f dp 
f 


/ (“g dp^) 


p « ^d 


and M^, 1 d^, and b^ as defined in Eq. (5.9) 


Momenttun balance for the natural circulation loop 
consisting of the steam drum, the downcomer and the riser 
yields the flow rate in the downcomer as 


w. = t20( P w| - 


w 


— -cvv ““ p "^ “p ^ 1.056+0.274 ■f (1+ 


2400 } + 0.95z] /5.51 


(5.11) 


where 


= 2V^ (b^ P + ^2 hj.) - 2V^ Pp a3_(fiy-a2P) P 


= Ip (^fgAfg) 

^ _ d , 

^2 ~ dp -^f 

2 / • 

and w^ * w^ +p ^ (b^ P + ^2 

where b^ and b^ are given in Eq. (5^9). 

Details of the formulation are provided in Appendix C. 

5 . 5*4 ^thermodynamic Properties of Steam 

Over the design range of pressure from 500,psia to 
700 psia, the specific volume of saturated liquid is assumed 
to be a constant whereas other thermodynamic properties are 
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assumed to vary linearly wit]! pressure. The functional 
relationships are provided in Appendix 0. 


3 •3*5 Steam Plow to Turbine 

Sonic flow through the valve is assumed as a convenient 
way of eliminating the consideration of downstream pressure 

34 

effects . The flow through the valve is then dependent 
upon the valve area and the upstream pressure only and the 
steam mass flow rate is written as 


Wg = constant .xa^ ( PAg )'^ 
where a_ = area of steam valve opening. 

3 

Normalising the flow with respect to the rated opening 
of the valve, the steam flow can he written as 


Za... = V' = a [ -.. P .-... ]* 

^oo ® ® Pioo - 


( 3 . 12 ) 


where a = normalised flow a.rea, w’ = normalised steam mass 

o 3 

flow rate, and P^qO^ ^g,100 pressure and specific 
volume of saturated vapour at hundred percent load. 

3-3*6 Pressure of Primary Heat Transport System 

The primary heat transport system is of constant 
volume with coolant in the liquid phase. The pressure of 
primary cooljknt in the outlet headers of the reactor is 
assumed to djepend on the mass of coolant in the entire 
circuit and/ the average temperature in the circuit . The 
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average temperature is calculated as a weigiited value of 

35 

the temperatures in the reactor and the heat exchanger"^ 

> 

and is given "by 

= 0.2126 + 0.2069 + 0.0862(12+^3) 

+ 0.2010 1^ (3.13) 

With Vp being the volume of the primary circuit, the pressiire 

35 

of coolant P is given by 

1 =C ) - C ] X 14.22 (3.14) 

and = ■V'd/Mjj 

where P = pressure in, psia 

= specific volume (in cm^/gm) 

Yp = volume of the coolant circuit 
= 4.28x10'^ cm^ 

Mjj = mass of the primary coolant in the circuit (in gms) 
A' = 0.28597 \ - 0.6909 

B» - 8594.06 - 111367.06/1^ + 365541.88/1^ 

C‘ = 16792.78 - 242883. 66/1^^ + 865909*5/^ 
and Ij^ = 0.01 [ ( lp^^-32)/1.8 + 273)] 

It is assumed that the mass of the coolant can be 

controlled by a feed or bleed valve such that 

as constant x a^ 


w 


(3.15) 
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where ajj = valve area normalised with respect to the 
total valve area 

In the above equation a^ is expressed as a percentage 
and Mjj in kilograms. I^e constant is taken as 0.117 which 
corresponds to a mass flow of 11,7 kg/sec (volume flow of 
817 litres per minute) at rated valve opening^^, 

Ihe feedwater flow to the inlet of tjie preheater leg 
is controlled by a valve, jt is assumed tiiat the flow through 
the valve is dependent on the valve area and is independent 
of the drum pressure. Therefore the flow rate is written as 

X area of valve opening, where is a 

constant. 

Normalising the flow with respect to full flow at the rated 
opening of the valve, the feedwater flow can be written as 

i '*^mi ^m ^ ■^m,100 

where v is the rated flow at full load 

m,iuu 

and a^ is the area of the valve opening expressed as a 
fERCtion of rated area at full load 

Pur):her,the transport delay between the reactor and 
the heat exchanger can be expressed as 

Tl(t) = '> 

and = 2 - 5^i(t) 


( 3 . 16 ) 
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Here refer to the temperatures of the coolant 

at inlet to and outlet of the reactor. and refer 
to the temperatures of the coolant at inlet and outlet of 
the heat exchanger. 

Equations (3»l) to (3.16) can he arranged in the 

form 

X = f(x(t), u(t), w(t), x(t~ Tj_), x(t- T 2 )) 
and y = g(x(t), u(t), w(t)) ( 3 . 17 ) 

where 

x^ = (n', c^,i=l, 2 , ...-, 6 , 

P» yr Mj)) 

J = (6k,, a^) 

and = (p, P, y) 

Here x refers to the variables which characterise the system, 
u refers to inputs which are controllable (or which are 
used for controlling the performance of system) 
w refers to inputs which are uncontrollable 
and .7 refers to certain outputs of the system • 

Equation 3*17 then represents the dynamical behaviour of the 
plant* 

3-4 STE4DY STATE GONhITIOHS 

m order to study the djmamic behaviour of the reactor; 
around full load conditions, steady state values of X and y 
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at full load conditions namely x^and y;^res|)ectively, are 
needed. These are obtained by the following procedure: 

If the values of u and w at full power operation of the plant 
namely u°and w° respectively, are specified, then Eq. (3.17) 
can be integrated till the steady state is reached. This is 
equivalent to obtaining the transient response of the plant 
for specified ]x and w from some initial x. In this study 
the initial value of x is formed from the values obtained from 
the steady state studies of the various sub-systems such as 
the riser, boiling and preheat legs of the heat exchanger and 
the natural circulation loop. The values of u°, w° and other 
parameters which characterise the full power operation of the 
reactor are given below. 


i) The reactor is operating at full power, i.e., in 

Eq. (3*3) is taken as 1.0. 

ii) Since AT„_ and AT „ at steady state are to be zero, 
6kfb in Eq. (3.5) is taken as zero, 

iii) 6k^ = 0 

iv) To maintain the water level at the centre of the dnun, 

feed water flow rate w^^^^ is taken to be equal to the 

steam flow rate w . Since the flow rates -of steam and 

s 

feed water are equal at full power, this leads to a_=a_. 

Hi jS 

v) a^j is used to adjust so that P is held at 1250 psia. 

vi) Peed water inlet temperature T is taken as 340®P. 

mi 
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vii) Steam valve opening is at the rated design 

value, i.e., a„=1.0. 

s 

With, these values of u° and w^, ^ and ^ are 
obtained as 

(x°)^ = (1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 

1075.6, 520.7, 503.4, 491.6, 550.65, 

492.3, 488.6, 592.3, 480.7, 416.9, 

474.7, 61.6, 568.9, 0.0, 0.363x10^) (3.19) 

= (568.9, 520.7, 1250.3, 0.0) 

Also u° and can be written as 

(li°)^ = (0.0, 1.0, 0.0) 

and (w°)^ = (1.0, 340.0) 

These values of and w^ represent the 

steady state conditions of the reactor at full power operatioii. 
Further these values give a steam flow rate of 8.714 lb /sec. 
per heat exchanger. This corresponcls to a total steam prod- 
uction of 2.5096x10^ Ib/hr. at 569 psia. 

The major advantage of obtaining the steady state 
values by a transient response analysis, as compared with 
procedures computing the steady state conditions by purely 
time-independent relationship such as minimisation of a 
function, is as follows There is a need for the verification 
of dynamic equations and this can be achieved by bringing the 
plant to steady state through a dynamic analysis. Obtaining 
a satisfactory steady state check is an indication of the 
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absence of gross errors in the model. 

3.5 TRANSIENT RESPONSES 

Effect of various inputs, both controllable and 

uncontrollable , on the system which is initially at steady 

ftill load operating conditions is now studied. A number of 

Inputs such as change in steam valve area a , change in the 

s 

inlet temperature of the feed water T oscillations in 

mi 

input reactivity 6 k , change in the area of feedwater valve 

c 

a are applied and the responses are obtained. However for 

in 

the sake of brevity, the responses due to two inputs namely 
change in a^ and change in T^^ are presented in this section. 

3.5.1 Response for a Step Change in Steam Valve Opening: 

Figure 3 . 1 (a-g) shows the response for a -5 percent 

# 

step change in valve area. 

The step change causes a sudden decrease (Fig. 3 *le) 

Ln the steam mass flow rate (w’ ), which initiates the trend 
of pressure rise in the drum (Pig. 3. Id). Increase in pressure 
leads to an increase in saturation temperature. This, with 
a constant temperature for the incoming primary coolant, 
results in decreased heat transfer and reduced steam production. 
Prarther, the reduced energy transfer from the primary coolant 
leads to a rise in the temperature of the primary coolant 
leaving the heat exchanger. However, this effect is felt 
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at the reactor only after a period equivalent to the transport 
delay* Reduced steam^-production coupled with an increase 
in steam outflow which is due to the pressure rise, reduces 
the rate of pressure rise with the pressure eventually attaining 
a steady value (Fig. 3. Id). 

Eie steam flow' is determined hy the flow area and the 
steam pressure. The sudden decrease in the steam mass flow 
rate, following the closure of valve, cannot he maintained 
as pressure increases. Hence the steam flow gradually rises 
with the press\ire, reaching a final steady value which is 
slightly more than the original equilibrium value (Pig,3..1e). 

The drum level is directly related to the mass of 
water in the drum which in turn depends on the flow rates 
entering and leaving the dr\am. With an increase in pressure, 
the quality of steam in the flow from boiling and riser legs 
decreases. This increases the amount of liquid held in the 
legs leading to a decrease in the mass flow rates of wat^r 
leaving them. This accounts for the initial drop in drum level 
(Fig.3*lg)« However, the decreased steam flow rate from drum, 
with a constant feed water flow to the preheat leg causes 
the drum level to rise in a short time. The rise continues 
for 130 seconds until the steam flow rate is more than the 
feed water flow rate, after which the drum level drops 
steadily* 
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Trxe increase in drum pressure leads to a rise in the 
temperature of the primary coolant entering the reactor. With 
the external reactivity input being held at equilibrium level, 
this increase leads to a rise in the average temperatxrre of 

the coolant (P^g.^olc) and a consequent rise in the temper- 

/ 

ature of the 'coolant leaving the reactor, ^his results in a 
small positive reactivity which raises the powerlevel of the 
reactor. However, the increase in fuel temperature due to the 
increase in temperature of the coolant and in the power level 
introduces negative reactivity to compensate the excess 
reactivity, £ind the reactor power reaches a steady value 
(Eig.3.1a). The average temperature of the primary coolant 
is higher than the equilibriiam steady state value. Since the 
mass of coolant is unaltered (the feed /bleed valves are not 
controlled), the primaiy circuit pressure rises enormously 
(Pig. 3 . If) before attaining a steady value. 

3 . 5.2 Response for a Step Change in Feed Water .Temperature: 

Figure 3«2 shows the response for a-15 ®F step change 
in the temperature of feed water entering the steam generator. 

The step change in feed water temperature reduces the 
average temperature in the preheat leg, leading to an increase 
in energy transfer to feed water and a consequent decrease 
in the temperature of the primary coolant leaving the steam 
generator* The increased energy transfer to feed water causes 
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the Gntha.lpy at the outlet of the preheat leg to remain 
constant or increase slightly. Hence immediately after 
the application of the step there are no changes in drum 
pressure and in drum level (Pig.5.2d and 3.2g). The drop 
in the outlet temperature of primary coolant at the steam 
generator is felt at the reactor after a time lag equivalent 
to the transport delay (Fig. 3 , ah). With the external 
reactivity input being held at its equilibrium level, the 
reactor power is constant at its equilibrium value. This, 
coupled with the reduced inlet temperature of the coolant 
to the reactor, results in a decrease of the outlet temperature 
of the coolant. Hence after a time, which is eq^ui valent to 
the Slim of transport delays between the reactor and the 
steam generator, the steam generator is supplied with the 
coolant at a lower temperature. This causes a reduction 
in steam production leading to a decrease in the drum pressure 
and the steam flow (Fig.3»2d and 3»2e). Further, a reduction 
in steam production increases the density of steamwater 
mixttare leaving the recirculating and boiling legs and the 
mass of liqtiid contained in these legs. This leads to a 
decrease in the flow of liquid to the steam drum. This in 
turn is reflected (Fig.3*2g) as a decrease in the liquid 
level. With the feed water mass flow at its equilibrium 
value, the reduced steam flow soon explains the rise in the 
liquid level. 
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xb-e decrease in the average temperature of the 
coolant results in a small negative reactivity which reduces 
the power level of the reactor. However the decrease in 
fuel temperature due to the decrease in the temperature of the 
coolant and in the power level introduces positive reactivity 
to compensate the negative reactivity, and the reactor power 
reaches a steady value (Fig, 3* 2a). The average temperature 
of the primaiy coolant is lower than the equilibrium steady 
state value. Since the mass of coolant is unaltered (the 
feed /bleed, valves are not controlled), the primary circuit 
pressure decreases enormously before reaching a steady value 
(Fig*3.2f), (Eq.3»14 is assumed to hold good, i.e., the 

heavy water is in liquid phase only), 

3.6 CONCniSIOHS 

These responses appear realistic and the trends 
compare with the expected behaviour of the plant. The 
response of the steam pressure and water level in the drum 
agree qualitatively with the results obtained in field tests 
of a boiler ' and with other theoretical results on drum 
type natural circulation boilers^* , However the proposed 
model can only be validated by comparison with the \incontrolled 
responses of the plant. Unfortunately such test data are not 
available. Since the nature of the responses agree with 
physical reasoning, there is added confid^ce in the model. ' 




49 


In tills chapter, a nonlinear model for the 
pressurised heavy water reactor is formulated. Uncontrolled 
responses of the system for various disturbances are obtained 
and are discussed. In order to design a control system for 
the reactor, the nonlinear model is linearised around oper- 
ating conditions. This results in a linear time invariant 
system which is valid for only small deviations from the 
operating point. This linear model is developed in the 
next chapter. 



CHAPTER 4 
LI HEAR MODEL 

4.1 INTRODUCTIOH 

Many of the differential equations developed in tlie 
previous chapter are nonlinear. Eor small perturbations 
about a steady operating point, these equations can be 
linearised to provide a set of linear first order ordinary 
differential equations. Such a system will facilitate the 
use of well developed control techniques in designing feed- 
back controls. Hence in this chapter, the plant is formulated 
as a linear system. 

4.2 MOLIPICATIOITS OP THE HOULIHEAR MOLEL 

The nonlinear model of the plant developed in Cha^pter 
3 is modified as given below. 

A simplification is made in the neutron kinetics. 

Instead of six groups of delayed neutrons, only one average 
precursor group will be used in subsequent analysis. The 
decay constant of the one group precursor is then calculated as 

, ^ _g 

^ 6 

i: I^i/^i] 

i=l - 

= 0.0708 sec"^ 

Transport delays between the reactor and the steam 
generator are represented by a first order approximation such tha 
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Since 


\ \ y -^2 


T = 2 T — T 
ro , ac ri» 


T. 


(2 ’^ac- - *■.)/ 


ri 


(4.1) 


U.2) 


Here lpj_» 35^^ refer to the temperature of the coolant 
at the inlet and at outlet headers of the reactor respectively; 
^1* ^4 temperature of the coolant at the inlet 

and at the outlet of the steam generator respectively. 


Ihese Eqs.(4.1) and (4.2) along wi^h the equations 
developed in the previous chapter are considered to represent 
the system and are linearised around the steady state values. 

4.3 IINBAHISATIOH PROCBDURE^'^ 


Considering the general equation of the form 

f(Xi, U2,...Uj^, x^, x^.-.x^) = 0 

and perturbing about the steady state values to eliminate the 
nonlinearities, gives 
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variables are evaluated at the initial steady state values about 
which the dynamic behaviour of the system is to be analysed. 

The unknowns in the equation are the pertiirbed values 

iiu^ > • • ^^1’ • • • ^ 

4.4 LIEEAR morel 


Eor example, considering the neutron kinetics . equation 
one average precursor. 


dJil ^ <5^ + -1 c' 

where 6k =. + “c(^ac-<c> (^^c'^^c> 

Initial steady state values are: T ~ 


af 


T 


af 


T = 1° 
ac ac 

6k = 6k° 
c c 


n' 


= 0.0 


= 1.0 
= 1.0 


Pert urbing the equation by the procedure given in Section 4. 3* it 
can be written as 


dt + (-^ )^c' + 

JO JO^ 

“f^^af ( ' * )o' ®c ^ \c * ^^c 

a ^ 

where suffix o indicates that the. terms are evaluated based 
on steady state values and regarded as constants in further 
analysis. Hence this can be written as 

A n* *= a^ An'+Eg + >5 '^^af ^ ^4 ^ ^ac ^5 *^^0 
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This is now a linear equation in the perturhed variables 

Ac’, 

In section 4»2, modifications are proposed to the 
nonlinear model developed in Chapter 3. These modifications 
are incorporated in the dynamics of the plant namely Eq.(3*17) 
and the resulting equations are linearised by the above 
procedure. A set of linear ordinary differential equations 
is thus obtained and this can be expressed as 

k = A S + 111 + ^2 - (4*3) 

where x, u, w represent, hereafter, the perturbed values i.e., 
AX, Au, AW respectively, from their steady state values, x, u 
and w are vectors of 19, 3 and 2 variables respectively and 

At I 2 matrices of orders (19x19), (19x3) and (19x2) 
respectively. Thus for the power plant, 

S® = (n,o, Ij, T^ 2 > % 3 ’ 

V' "p’ "d’ P’ 3^’ “d) 

and H* = (ag, \i) 

where n = changes in the power level expressed as a percentage 
of full power 

y = changes in water level expressed in the unit of 
0.01 inch 

6k = reactivity inputs in cents, with 100 cents being 
equal to a reactivity of 0*0064. 
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^m* changes in the ariSa of the feed water control 

valve, primary circuit pressure control valve area and steam 
valve area respectively and are expressed as percentages of 
full opening (positive and liegative values for a^j refer to 
the feed and the bleed val^-es respectively) and the other 
symbols having the same meaning as in Chapter 

Further, Eqs. (3»14) and (3.12) for the calculation of 
the primary circuit pressure and the steam mass flow rate can 
be linearised in a similar manner and this yields 

P = 2.214 + 49.14 - 0.64 x^g + 23.93 x^.^ 

+ 9.963 x^+ 9.963 Xg + 23.24 x^^ (4.4) 

and 

w^ = 0.1776 x ^5 + (4,5) 

where x^, u^, w^ refer to the i"^^ variable of the vector 
X, u and w respectively. 

Prom Eq. (4.4), 

P = 2.214 + 49.14 x^.- 0.64 x^ + 23-93 x^^ 

+ 9.963 X 5 + 9.963 Xq + 23-24 

Substituting for the derivatives of the state variables x^ 
iErom Eq. ( 4 . 3 ), 

P =10.8 x^ - 156.0 x^ - 25.9 x^ + 27.7 Xg + 23.7 Xg 
+ 23.1 ”” 67.0 Xj^j^ 9.89 ^^2 5.33 

+ 151.0 x^Qt 0.262 u^ 


(4.6) 
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Bq. (4«6) is then used in place of the first order differential 
equation for in Eq. (4*5) und Eq.(4»5)' can "be written as 

X = Ax + B^u + B^w ( 4 . 7 ) 

with the last variable x^g now being the primary circuit 
pressure P, 

The matrices A, ^ and B 2 for the above system are 
provided in Tables 4.1, 4.2 and 4 . 3 . 

4.5 TRANSIENT RESPOITflBS 

Uncontrolled dynamic behaviour of the system represented 
by Eqs.(4.7) and (4-5) is obtained for two inputs, which are 
same as those considered in Chapter 3, namely a -5 percent step 
in the area of the steam valve w^ and a -15 degree' step in the 
feed water temperature w^, with the initial conditions being 
x=0 and u=0 at time t=0. These responses are shown in Fig.4»l 
and 4 . 2 . These transient responses are in close agreement with 
those obtained for the nonlinear model (in Chapter 3), thereby 
validating the linearisation prcdedirre. The difference between 
the responses shown in Fig. 4 . 2b, c and Fig. 3.2b, c is mainly due 
to the first order approximation of the transport delays 
between the reactor and the steam generator. A higher order 
approxination will improve the agreement during the early 
part of the transients. But. this will increase the oomput* 
'ational efforts needed for control studies. In view of this 
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Table 4.2 

Matrix (19x3 Matrix) 

^ has 7 non-zero elements and 50 zero elements. 
The non-zero elements are; 


1, 

B^(10, 

Bi(12, 

Bi(13, 

\(15, 

Bl(16, 

B^(19, 


1) = 0,942+01 

2) = -0.17E+00 

2 ) = - 0 . 332+00 
2 ) = -0.422-01 
2 ) = 0.732-02 

2) = O. 3 OE+OO 

3) = 0.26E+00 


Table 4.3 


Matrix B2 (19x2 Matrix) 

B^ has 4 non-zero elements and 34 zero elements. 

The non-zero elements are: 

B2(15, 1) = -0,18E+00 
B^tie; 1) = 0.552+00 

B 2 ( 12 , 2 ) = 0.292+01 

22(15, 2) = -0.13E+G1 

In the tables given above A(I,J) represents a. .of 

I ' ' ' ■ . , , X- f ■ 3 ' ■ ' ' ' 

matrix A* Also, the symbol E stands for the power of 

10. Ecr example A(l,l) = -0.942+01 indieates ^=0,94x10'*'^. 
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tiie first order approximation is considered as acceptable. 
Hence the system represented by Eq*(4.7) can be utilised 
for studying the dynamical behaviour for small perturbations. 
In the sub frequent chapters viz. Chapters 5 and 6, the 
linear system is studied to develop feed back controls based 
on the optimal control theory. 
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CHAPTER 5 
OPTIMAL CONTROL 

5.1 INTRODUCTION 

The previous chapter dealt -with the \mcontr oiled 
hehaviour of the plant when it is acted upon hy disturbances* 

It can be seen that the effect of the disturbance was to vary 
most of the system variables from their initial values* Ihe 
basic requirement of the power station is to supply electrical 
energy as demanded by the power distribution network; this 
necessitates the control of a few important plant variables 
like the temperature of the coolant in the reactor* the 
pressure of the primary coolant circmt and the water level 
in the steam drt© at prescribed values,. 

The control of the Variables in the power station is 

12 6 

presently accomplished by analog devices ’ which are usually 
designed to respond to an error in a single variable. A 
number of independent controllers are used to regulate variables 
such as water level in the boiler, steam pressure etc. The 
design of such standard controllers neglects the effects of 
dynamic interaction existing between the controlled variables^ 
that is, if an error occurs in one of these variables, its 
controller moves a single actuator. Unfortunately, movemoit 
of this actuator causes errors in other controlled variables 
due to interaction among them. The new steady state operating 



66 


point is reached after a number of controllers haisre acted 
in an independent manner. 

In alternative approach to plant control is to design 
a single multivariable controller; input to this controller 
consists of the desired and the actual values of the variables. 

If an error occurs in a plant variable, the controller responds 
by moving a number of actuators in a co-ordinated manner so 
as to reduce the error without causing \uinecessary errors in 
other variables and thereby enabling the system to reach 
quickly the now steady state. The application of digital computer 
in power station operation for data logging, for alarm scanning 
and for normal start-up and shut-down, introduces the possibility 
of developing and implementing a multivariable controller for 
the power station. In this chapter, optimal control theory is 
applied to the linear model of the pressurised heavy water 
reactor to design a multivariable controller. 

5.2 BRIEF REVIEW OF OFTIMAI COETROL THEORY 

Analysis and design of linear multivariable controller 
by using optimal control theory have been well discussed in the 
literature^®* The objective of the multivariable controller 
is to generate a linear constant state feedback control law 
which will transfer the system from a given initial state 
to the desired final state. In so doing the control system 
must satisfy the requirements relating to the performance oi^^ ^ 
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the system such as the desired output, desired control effort 
etc. , and also the implementation. The requirements on the 
performance is represented mathematically as a performance 
index. Quadratic performance index in the output or state 
variables and the control variables ■which minimises the error 
in the output or state variables and the control effort is 
most frequently used. Tor a linear system, the quadratic 
performance index yields a linear control law which is 'easy 
t O CoYYvjbU/te , 

The performance criterion is hence taken as 

j i £(t^) + i [x^(t)a x(t) + U^(t)gu(t)] dt 

( 5 . 1 ) 

where x and u refer to state and control vectors,. 
jQ is a positive semi-definite matrix, 

R is a positive definite matrix, 
t^ is the initial time, 
and t^ is the final time • 

A unique set of weighting matrices, i.e,., the matrices 
R and S to satisfy the prescribed or required design 
specifications generally does not exist, However, the lack 
of uniqueness of the weighting factors does introduce a 
flexibility which makes selection of the performance index 
simpler and meaningf\il. Depending upon the relative importance 
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of the state or output variables and the control variables, 
the cost matrices jQ and R are chosen to reflect the desired 
closed loop performance of the system. The constraints on 
the control and output variables can be indirectly taken into 
account by assigning suitable penalties to the constrained 
variablespand ^ and R matrices can be chosen to reflect 
these pensilties. For example, a large value for ith diagonal 
element of R means that the corresponding control variable, 
i.e. ith control variable, is not allowed to vary significantly. 

Having selected the performance index, the design of 
multivariable controller can be stated as: To find the control 
u which minimises the performance index subject to the constraints 

;^ = Ax+Bu 

and X ("tQ) = Sq 

where x is a (nxl) vector of state variables, u is a (mxl) 
vector of control variables, and A and B are real, constant 
matrices of appropriate sizes. The control which satisfies 
the. above requirements can be obtained by the application of 

38 

Pontryagin maximum principle or the Hamilton Jacobi equation 
and is well known in optimal control theory. The existence of 
the optiml control however, requires the prior investigation 
of the system for controllabilit3r^^ as discussed in Appendix 
D. The optimal control can then be obtained (as shown in 

1 IT "R'i n<={ 
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li(t) = £(t)x(t) (5.2) 

where P(t) is the positive definite solution of 

i (t)= -p(t) A - A%(t) + P(t)B B^ t(t) - a (5.5) 
with a terminal condition given hy P(tf) = S* , 

The above equation, i.e., Eq.(5.3) is known as the matrix 
Riccati equation and the solution of this equation provides 
the matrix P to calculate the feedback control law through 
Bq .(5.2). This problem of finding optimal control policy 
for a finite t^ is known as the finite time linear regulator 
problem. 

When t^ in Eq. (5*1) is infinite, P can be obtained by 
solving tho algebraic matrix Riccati equation, i.e. the 
solution of • . ■ 

P A + - P B R~^ B^ P + ^ = 0 

gives the matrix £ to calculate the control law. Also P is 
independent of time. This problem is known as infinite time 
linear regulator problem. The procedures for solving the 
algebraic and the differential matrix Riccati equations are 
discussed in the following sections. 

5.2.1 Infinite Time Regulator Problem 

Infinite time re giiL at or problem requires the solution 
of the algebraic matrix Riccati equation which is obtained 
by the method of successive approximation^*^^^^. In this 
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method, an initial choice of P is made such that the controlled 
system given hy the equation 

i = U - S P) X (5.4> 

is stable. This initial value of P nay he guessed or may 
he obtained by integrating Bq. (5*3) backwards in time till a 
P which makes the controlled system stable is obtained. 

Then an iterative scheme is followed to obtain P at the ith 
iteration from its value at the (i--l)th iteration. At the 
ith iteration, P is obtained as the solution of 

+ ( 4 ^)^ = 0 

where 4^ = 4 ~ 1 P^“^ 

i i— "I —1 T i— 1 

and S = ^ £ B R B^ p-^ 

These iterations are carried out until the difference between 
the elements of successive P matrices is within the required 
accuracy. In this work, the iterations are continued till the 
elements of P matrices at successive iterations agreed upto 

5 

1 part in 10 , 

5.2.2 Binit e Time Regulator Problem 

This requires the solution of the matrix Riccati 
equation in its differential form, i.e. , Eq. (5.3). The. 
solution is obtained by direct integration of the equation 
backward in time, starting from the final time t^ with B(t^)=S 
till the! initial time t^ is reached. The matrices P(t) thus 
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obtained at various time intervals, are stored and are used 
for generating the feedhack control given by Eq.(5.2). It 
can he seen that the finite time regulator problem results in 
a time varying control law; the implementation of the feedba.ck 
law then requires a programmed control and becomes difficult. 
Therefore in this thesis, this is not attempted and the 
infinite time regulator theory is used to obtain a time 
invariant control law for the pressurised heavy water realtor. 

To sucimarise, the problem of calculating the feedback 
control law for the pressurised water reactor can be treated 
as a linear regulator problem. The infinite time regulator 
problem which yields a time invariant control law, can be 
stated as; Given 

X = S + Ji f 2 ( 0 ) = Xq 
it is required to find the control u which minimises the 
performance index 

J = x(t) + u^(t) R u(t)] dt 

where ^ and R are appropriate matrices. The solution to the 
above problem is written as 

u(t) = K x(t) 

= 1 x(t) 

-IT 

where E is the state feedba^ik matrix and is equal to -R £ R 
and where is the positive definite symmetric matrix 
satisfying the following matrix Ric cat i equation 
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£ii + A^Il-£Bgr^£,^P+^ = 0 , 

Tlie state feedback matrix g is then used to generate the 
control vector u(t) from the state vector x(t). 

5.3 CONTROL STUDIES ON THE PRESSURISED HEAVY WATER REACTOR 
The infinite time regulator theory described in the 
proceeding section is applied to develop proportional controller 
and proportional-integral controller for the pressurised heavy 
water reactor. The following sections describe the performance 
of the reactor with these controllers. 


5 . 3»1 Proportional Controller 

The dynamics of the linear model of the pressurised 
heavy water reactor has been represented by Eq.( 4 - 7 ) which is 
repeated below for convenience. 


X 


__ 4 ^ + %]i + S2- ( 4 * 7 ) 

The design and operational requirements of the reactor impose 
minimum variations of the following three state variables; 
i) average temperature of coolant ii) level of water in the 
steam drum Xj_g and iii) pressure of the primary heat transport 
system In order to meet these specifications, the 

weighting matrices ^ and R are chosen as 

^ = diag (0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 

t 0, 0, 0, 0.1,: 0, 0, ^ ^ ^ ( 5 . 5 ) 
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R = diag (1, 1, 1) 

where the perfoma,nce index J is given by 


(5.6) 


J = i J Ix^Ct) S£(t) + u^(t) R^Bj(t)3 dt 

The aininisa-tion of J using infinite time regulator, theory 
leads to 


u(t) = g x(t) 
where the elements of 
Table 5.1. 

Matrix 

i;(l,l) = -C.9973E-01 
E(l,2) = -0.9234E-C0 
K(l,3) = -0.8569E-00 
K(l,4) = +0.1164E+02 
Z(l,5) = +C.1237E+01 
K(l,6) = -0.3628E+00 
K(l,7) = +0.2222E+01 
E(l,8) = -0.3200E+00 
K(l,9) = -0.1356E+01 
K(l,10)= +0.2575E+00 
Z(l,ll)= +0.5238E+01 
K(l,12)= -0.1686E+00 
K(l,15)= -0.1719E+00 
K(l,14)= +C.5599E-01 
K(l,15)= -C.82C7E+00 
E(1,16)= -0.9063E-01 
K(l,17)= -0.1287E+01 
K(l,18)= -0.1142E+02 
K(1,19)= -0.3013B+00 


the feedback matrix g 

Table 5.1 

g : (3x19 matrix) 

K(2,l) = -0.1388E-02 
K(2,2) = -0.1713E-01 
K(2,3) = -0.9041E-02 
Z(2,4) = -0.4163E+01 
K(2,5) = -0.1530B+01 
K(2,6) = -C.2839E+00 
E(2,7) = +C.7054E+01 
E(2,8) = -0.1764E+01 
E(2,9) = ~0.9120E+00 
E(2,10)- +0.1432E+01 
E(2,ll)= -0.1852E+01 
E(2,12)= -0.1827E+00 
E(2,13)= ~0-4355E+00 
E(2,14)= +0.2657E+0C 
E(2,15)= +0.43G5E+0G 
E(2,16)= -G.3G82E+00 
E(2,17)= -G.3317E+C1 
E(2,18)= +0.1G27E+G1 
E(2,19)= +G,8652E-G1 


( 5 - 7 ) 

given in 


E(3,1) = -0.8248E-G2 
E(3,2) = -0.7908E-01 
E(3,3) = -C.7466E-01 

E(5»4) = +C.19vCl+Cl 

E(3,5) = +0.1196E+OC 
E(3,6) = -C.2G01E+GC 
E(3f7) = -C-5945E4-00 
K(3,8) ,.= +0.228CE+C0 
E(3,9) = -0.2917E-G1 
E(3,1G)= -G.1258E+G0 
E(3,ll)= +0.8560B+G0 
E(3,12)= +G.9244E-G2 
E(3,13)= +0.5C22E~01 
E(3,14)= +0.3455E-.02 
E(3,15)= -0.1099B+0C 
E ( 3 , 1 6 ) = +0 . 2381E- 01 
E( 3 , 17 ) = +0. 4150E+00 
E(3,18)= -0.1220E+C1 
E(3,19)= -0.4956E-01 


[The notation used to represent the elements of g is 
the same as that used in Chapter 4]. 
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I-erfornance of the system with this feedback control 

is obtained for a constant disturbance of -5 percent change 

in the area, of steam valve Implementation of this control 

is seen to produce stable closed loop responses (Fig. 5*1) but 

the system settles with undesirable stea.dy state off— sets. 

These off-sets depend significantly on the weighting factors 

used, Fuirther the off-sets are to be expected since the 

performance index does not take into account the physically 

necessary off-sets which must exist in the variables following 

10 

a sustained change in the disturbance vector . However the 
steady state off-sets can be eliminated by the method discussed 
in the following section. 

5.3*2 Proportional Controller Based on Error Coordinates 

Ihe steady state off-sets obtained in the previous 
section can be eliminated by modifying the performance index. 

A new performance index is formulated^ ^ as 

J = i ^ I (x (t) -Xg)^ S (X(t )-Xg) + (U(t)-Ug)^ 

S (]i(t) - lig)3 clt (5.8) 

In Ba.(5. 8) x and u are the equilibrium values after a 
^ ^ ' — S, —S,' 

step nhange in the disturbance vector. 

5. 3. 2.1 Calculation of Steady State Values 

The steady state values x^, u^ are determined from 
Eq.(4.7). i.e., the steady state of the system after a step 












change in w- can he obtained from Eq.(4.7) as the 
solution of 
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io = ii ^ + In + Bo w = 0 (5.9) 

Having specified w, it can be seen that the above equation 
has more unknowns than the number of equations. This allows 
for certain unimowns to be specified arbitrarily as required 
by the needs and design considerations. In this pressurised 
heavy water reactor study, there are 22 variables (19 state and 
5 control variables) in 19 equations. Hence three variables, 
either state or control, can be specified. Itirthermore, 
since there are three control variables, , three of the states 
can be driven to zero, i.e. to their initial values. These 
three states are chosen as the average temperature of coolant 
in reactor the level of water in steam drum x^g and the 
pressure of the primary heat transport system x^g and it is 
required to hold their magnitudes at their respective equilibrium 
levels prior to the application of the disturbance. Thus the 
steady state values of x^, x^g and x^g are specified as zero 
and Bq.(5.9) is solved for other unknowns. The solution then 


provides the necessary steady state values of the control 
variables so as to maintain x^, x^^ and x^g at their prescribed 

■ ' . 'jt 

values in the presence of a disturbance. 


For a —5 percent step in the area of steam valve w^, 
the above procedure yields the following values for x^ and ^ 
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=“(3.5196, 3-5196, 19.5617, C.O, -I.0310, -1.5318, 
3.0295, “1.4617, “1.6043, “5.6431, “1.4076, 

“1.1110, “6.3133, “1.0135, “9.5081, 0.0, 

1.4076, “1.4076, 0.0) 

=-(1.6815, 3.3011, 0.0) 

5. 3 .2. 2 Dynamics Based on Error Yariables 

The dynamics of the reactor is now based on the errors 
between the steady state values x and u_ and the Instantaneous 
values of x(t) and u(t). Equations (4.7) and (5.9) are combined 
to yield 

x(t) - X = A [x(t)“ X ] + B [u(t) “ u ] 

vD fc) Q 

Define vectors e and f as 
e(t) = 2 (t) “ Xg 
and f(t) = u(t) “ u 

Then the above equation can be written as 
^(t) = A £(t) + B^ f (t) 

and the performance index Eq. (5 . 8) is written as 

J = i ite^(t) ^ e(t) + f^(t) B f(t) ] dt. 

The application of linear regulator theory to the above system 
of equations with the same\Q and R matrices as given by 
Eqs. (5.5) and (5.6) leads to the following control law 

f = E e 
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The elements of the feedback matrix E are the same as that 
of Eq.(5.7) which are given in Table 5»1» Transforming the 
control law 

f(t) = K e(t) 

to the original co-ordinates x and u, the control u is 
writtifi as . 

liC't) = lig + S (s('t)-Xg) (5. 1C) 

The performance of the system with the above control law for 
a -5 percent step in w^ is shown in Pigure 5.2. It can be 
seen that this control law has produced closed loop responses 
with complete elimination of the steady state off-sets- in the 
selected variables x^, x^g and Ike response of x^, i.e. 

the average temperature of the coolant in the reactor, is 
obtained for a large period of time and. is compared with the 
response obtained with the control law given by Eq.(5.7) 

(Pigure 5.5)* Comparison of figures 5*1 and 5»2 show that the 
initial values for the controls, A and A a^j are quite 
different fnr the two control schemes, i.e. Eqs.(5.7) and 
(5.10), with the control law given by Eq.(5^10) requiring 
large initial values for Aa^ and Anj^. 

5. 5. 2. 5 Drawbacks . of the Controller 

The control law represented by Eq, (5.10) is theoretically 

■foA Mtfi 

optimal^ . However, it is quite difficult to iB^jlement because 
of the following two drawbacks. 
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(a) Equation (5.10) needs a knowledge of the nature of the 
disturbance and its magnitude so that the steady state vectors 

X and u can be evaluated. 

—s “S 

(b) This control law is sensitive to the parameters of the 
system. This is mainly due to the use of the terms (x - gg ) 
and (u - Ug) in the performance index which forces u and x 
to approach the calculated steady state values u„ and 
Since x and u are obtaitt.ed from Eq,(4,7)» their magnitudes 
depend on the matrices A, and B 2 » If the plant parameters 
differ from the nominal values, i.e., the matrices ^ 3-^ 

B 2 are different, then the control law given by Eq. (5.10) will 
drive the system to a steady state where the requirement of 
zero steady state off-sets in the selected variables may not 
be fulfilled. 

The shortcomings of the above control scheme are mainly 

due to the use of proportional controller in the presence of 

a constant disturbance. It can be visualised that a proportional 

integral controller may be able to meet the requirements of 

zero steady state off-sets in selected variables when the 

system is subjected to constant disturbances. Since linear 

regulator theory leads to proportional controller using state 

variable feedback, modifications of the regulator theory is 

42 

needed to yield a proportional integral controller . The 
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following section discusses the design of the proportional — 
integral controller using linear regulator theory. 

5.5*3 Proportional-Integral Controller 

In section 5.3*2, optimal control laws are calculated 

i . 

based on performance index which penalise the deviation of the 
state and control vectors from their steady state values. 

It is also seen that the steady state values depend on the 
nominal values of the plant parameters. If the plant parameters 
are different from the nominal values, then x and u would not 
converge to their calculated steady state values. This is a 
disadvantage with the proportional controller of section 
5»3.2. However this drawback may be eliminated by noting 
that X and u do tend to some constant values so that their 
time derivatives tend to zero regardless of the actual values. 
The way that this can be incorporated into the optimal regulator 
theory has been discussed by M. Athans^^, Parker'^^ and Smith 
and Davis on‘^‘^ and the procedure is now described. 

5. 3*3*1 Modified Linear Regulator Theory^^*^'^ 

Consider the system represented by 

S = + (5*11) 

i = gc + D^ u + D^ w (5*12) 

with x(o) = x^ , and u(o) ?= u^ (5.13) 

X is a (nxl) vector of state variables 
u is a (mxl) vector of control variables 


where 
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and w is a (qxl) vector of diBtur'bance.'varia'bles 

Further ^ is a (nncl) vector of output variables whose values 
are to be held zero. 4 , C, and B 2 real 

constant matrices of appropriate dimensions. 


Define vectors z and v as 






(5-14) 

(5.15) 


By differentiation, Eqs.(5.1l) and (5.12) are transformed as 


z 


= A z + B V 


(5.16) 


2 ( 0 ) 


where K = 


and B 


A Bjj^ B2 


C £1 £2 


A 0 
0 P 

D. 



Sc' 

i 

1^0 

i 

j 

w 


J 


(5.17) 



With the performance index defined as 

J = ■^r O ^ S E z] it (5.18) 

where Q and R matrices are of appropriate dimensions, a stable 
linear state feedback control law may now be obtained for 
the transformed system represented by Eq, (5.16). (The control 
is then of the form 
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s: = (5.19) 

which on transformation into the original variables and ■by- 
partitioning of the natPiJc ^ becomes 


k = E i + i Z (5.20) 

On integrating the above equation, one obtains 

t 

u(t) = E x(t) + 1 X ■ y: ds (5.21) 

o 

In obtaining Eq.(5.21), the initial conditions of x, u and 
y; have been taken as zero. Equation (5.21) provides a control 
law with proportional feedback of the state variables and 
integral feedback of the output variables which are to be 
driven to zero as time tends to infinity. The above technique 
is now used to obtain a proportional - integral controller 
for the pressurised heavy water reactor and is discussed in 
the following section. 


5. 5. 3. 2 Application on the Pressurised Heavy Water Reactor 

The dynamics • of the linear model of the reactor has been 

represented by Eq.(4.7) which is repeated for convenience 

♦ ■ ■ 

X = 4 S + ^ 21 + ^2 (4 •7) 

Since the design and operational requirements of the reactor 
need that the variations in the average temperature of the 
coolant in the reactor, the level of water in the steam dirum 
and the pressure of primary heat transport circuit are to be 
minimum, the output vector y becomes 

'' 

X ~ ^^4^ ^16^ ^19^ (5.22) 



87 


Hence for this "tlie matrices C, D^, in Eq.(5*12) are 
•such that all the elements are zero excepting 

°1,4 ^2,16 ®3rl9 ^ (5.23) 

The weighting matrices in the performance index given by 
Eq. (5.18) are chosen as 

^ = diag ( 0, 0, 0, Of 0, 0, 0, 0, 0, 0, 0, 0, 0, 

0., 0, 0, 0, 0, 0, 1, 0.1, 0.1) (5.24) 

and R = diag (l, 1, 1) (5.25) 


The modified linear regulator theory is now applied 

to the system represented by the equations (4.7), (5 . 23) » (5. 24) 

and (5.25). The control law thus obtained, is given as 

t 

u(t) = K x(t) + L J X ds (5.26) 

0 

where the matrices Z and 1 are shown in Tables 5.2 and 5.3. 

The above feedback control law is used to study the 


closed loop dynamic behaviour of the pressurised heavy water 
reactor for a -5 percent step change in the area of steam valve 


w^. The performance of the system is quite satisfactory (Eig,5.4). 
The output variables x^g and reach zero magnitudes 
within a minute of the occurrence of the disturbance. 


In order to study the effect of the weight on the 
average temperature of the coolant in the reactor, the matrix 
^ given Ifey Eq.(5.24) is modified as • 
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Table 5.2 

Matrix K: (3x19 Matrix 


K(l, 1) = -0.154E+00 
Z(l, 2) *= -0.990E+00 
Z(l, 3) = --0,139E+01 
E(l, 4) = +0.175E+02 
K(l, 5) = -0.308E+01 
E(l, 6) .= -0.458E+01 
K(l, 7) = +0.487E+01 
K(l, 8) = -0.158B+01 
Z(l, 9) = -0.297E+01 
E:(1,10) = +0.397E+00 
Z(l,ll) = +0.790E+01 
E(1,12) = -0,282E+00 
E(l,13) = -0.235E+00 
Z(l,14) = +0.770E+00 
Z(l,15) = +0.332E+00 
K(l,16) = -0.266E+00 
Z(l,17) = -0.421E+01 
Z(l,18) = -0.188B+02 
K(l,19) = -0.517E+00 


K(2, 1) = -0.426E-02 
K(2, 2) = -0.351B-02 
Z(2, 3) = -0.416E-01 
K(2, 4) == -■0.112E+02 
K(2, 5) = -0.180E+02 
E(2, 6) = -0.114E+02 
E(2, 7) = +0.267E+02 
E(2, 8) = -0.795E+01 
E(2, 9) = --0.629E+01 
E(2,10) = +0.342E+00 
E(2,ll) = -0.420E+01 
E(2,12) = -0.622E+00 
E(2,13) = -0.158E+01 
E(2,14) = +0.274S+C^1 
E(2,15) = +0.698B+01 
E(2,16) = -0.153E+01 
E(2,17) = -0.192E+02 
E(2,18) = +0.492E+01 
E(2,19) = +0.222E+C0 


Table 5.3 

Matrix I; (3x3 Matrix) 


1(1, 1) ■= -0.165E-I-00 
L(l, 2) = -0.811E-01 
1(1, 3) = --0.301E+00 


L(2, 1) = +0.794E-01 
1(2, 2) = -0.305E+00 
1(2, 3) = +0.778E-01 


h:(3, 1) 
e;(3, 2) 
h:(3, 3) 
E|[3, 4) 
£(3, 5) 
E(3, 6) 
E(3, 7) 
K(3, 8) 
E(3, 9) 
E(3,10) 
E(3,ll) 
K(3,12) 
E(3,13) 
E(3,14) 
K(3,15) 
E(3,16) 
E(3,17) 
E(3,18) 
E(3,19) 


1(3, 1) 
I-(3, 2) 
1(3, 3) 


= -0.142E-01 
= -0.352E-01 
= ~C.134E+00 
= +0,164E+02 
= +0.332E+01 
= -0.476E-01 
= ~0.233B+01 
= +0.326E+01 
= +0.259E+00 
= -0.293E+00 
= +0.761E+01 
= +0.331E-01 
— +C *1401+00 
= -O.llOE+00 
= -0.857E+0C 
= +0.103E+00 
= +0.653E+01 
= ~0*995E+00 
= -0.353E+00 


= 0.981E+00 
=+0.110E-Gl 
=-0.570E-01 


[The sotaition used to represent the elements of E is the 
same as that used in Chapter 4] 
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g = diag (a, 0, 0, 0, 0, G, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 100, 0.1, 0.1) (5.27) 
and the matrix ^ is maintained as in Eq.(5,25),, i.e. 

R = diag (1,1,1) (5.28) 

Optimal control law of the form given hy Eq. (5.26) is 
obtained with these weighting matrices and the response of 
the system for a -^5 percent step in w^^ is shown (as dotted 
lines) in Eigure 5.4. It can be seen that the increase in the 
weight on coolant temperatirre has reduced the variations in 
AT considerably. The reactivity 6k has changed marginally 
when compared with the reactivity needed for the previous 
weighting matrices. The magnitude of Aa^j has increased. This 
is due to the high sensitivity of the pressure of the primary 
heat transport system to the variations of temperature in 
the circuit. 

■» , 

5.4 CORCIUSIOES 

In this chapter, optimal control theory has been 
applied to the dynamics of the linear model of the pressurised 
heavy water reactor to obtain feedback controllers. The 
application of linear regulator theory has resulted in 
proportional controller and the closed loop behaviour of the 
reactor with the proportional controller is stable with steady 
state off-sets in the presence of constant disturbance. These 
off-sets are eliminated by modifying the performance index. 
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The modified performance index also yields a proportional 
controller hut it is not feasible to implement this controller 
since it requires a knowledge of the natiare and the magnitude 

Attf> 

Of the ^disturbance. 

In order to overcome the drawback of the proportional 
controller, a proportional -integral controller is envisaged 
and the linear regulator theory is modified to yield the 
proportional- integral controller and the performance of the 
system with this controller is obtained for a constant distur- 
bance. Eventhough the response of the system is acceptable, 
the implementation of the controller requires a knowledge of 
all the states. However only a few states like the average 
temperature of the coolant in the reactor, the water level 
in the steam drum and the pressure of the primary heat transport 
system are directly measurable and have practical significance. 
Hence it is pertinent to have a control law which is not necess- 
arily optimal, but is dependent only on the available states. 

A control of this nature, known as sub-optimal control, is 
developed in the next chapter. 



CHAPTER 6 


SUB- OPTIMAL CONTROL 

6j.1 IHTRODUCTIOR . 

The concept of state space and the use of Pontryagin 

maximum principle have resulted, as shown in the previous 

chapter, in optimal feed hack control la^i^^s which call for a 

complete measurement of the state vector. In many physical 

systems, the state variables need not necessarily correspond 

to physically measurable quantities. To realise the control 

law, it is often suggested that the unavailable state variables 

can be constructed via a ifelman-Bucy filter or an observer. 

But it is quite impractical to reconstruct the state variables 

in a large multi-variable ■ system such as those which occur in 

chemical processes, power systems etc. In these cases it 

would be convenient if the designer has the freedom to select 

based on experience or intutiop certain variables and make 

the control a function of these state variables* This may 

result in the deterioration of the performance index; but if 

the performance of the system is close to the performance of 

the optimally controlled system, then it would be advantageous 

to go in for this control law. This type of feedback control 

scheme where the control variables are linear combinations of 

45 

selected state variables is obtained in this chapter 
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The design of the constrained controller then reduces to the 
selection of values for the non-zero elements of the feedhach 
matrix E. 

6.3 OPTIMAL DESIGN OE CONSTRAINED CONTROLLER 

If the constrained feedback matrix K is to generate 
optimal control of the form u = K x for the system represented 
by Bq.(6.2) and the performance index J given by Eq.(6.1), then 
it is necessary that 3J / 9 kj^^ must vanish for each allowed 
feedback kj_j , Define a matrix G of order mxn with its element 
g. . representing the gradient of J with respect to the element 

ij Qk . . 

Then for optimality, g^^ must be zero for each allowable 

45 

feedback i and 3 . Jameson has shorn that the gradient matrix 


G can be written as 

G = 2 (B^P + R I) ¥ (6.4) 

where ¥ and P satisfy the following equations 

I I + I + So So = (^• 5 ) 

£ + £ £ + ^ 1 I = 0 (6.6) 

with £ = A + B E (6.7) 
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Further the value of the performance index J is given hy 

J = x^ P X 

—0 o 

= trace [P (6.8) 

Thus Eqs. (6,8) and (6,4) provide explicit expressions for 
the magnitude of the performance index and its gradients 
with respect to the feedback matrix K, The design of optimal 
constrained controller now reduces to the selection of the 
allowable feedback elements k. . such that g- .=0 for each 
allowable i and j . 

6 . 3.1 Elimination of Dependency on Initial Conditions 

It can be seen from Eqs. (6.5) and (6.8) that the 
gradient matrix G and the value of the performance index J ■ 
depend on the . initial condition x^, specified by Eq.(6.2). To 
overcome the dependency on initial condition, it is generally 
assumed^ that x^ lies on the surface of a unit sphere in the 
n dimensional eucledian space and x^ x^ is taken as an identity 
matrix I. Substituting I for x^ x^ in Eqs. (6.5) and (6.8) 
the matrix G and J can be obtained from the following equations. 

a) G is given by 

g = 2 (B^P + R K) ¥ (6.9) 

where ¥ and P satisfy the following equations 

P ¥ + ¥ + i = 0 

Z + E I + a + R K = 


0 


( 6 . 10 ) 

( 6 . 11 ) 
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■with 1 = . 4 + B E (6.12) 

Id) J is given by 

J = trace [P] (6.13) 

Procedure for obtaining the optimal feedback matrix 
E using the above equations is discussed in the following 
section. 

6.3»2 Method of Solution 

Having obtained explicit expressions for the 

performance index J and its gradient with respect to the 

feedback matrix K, the procedjure of obtaining the optimal 

controller is stated as: Selecting an initial guess for Z, 

Eqs. (6.10), (6.11) and (6.12) are solved^^ to obtain the 

matrices P and ¥. These matrices are then used in Eqs. (6.13) 

and (6.9) to calculate the performance index J and its 

gradient matrix G. Since the problem of designing the optimal 

controller is one of obtaining the elements k. . for each 

ID 

allowable i and d » which will result in the corresponding 


= 0, the problem is treated as one of minimising a function 
using its gradients. The function thus becomes the perfor- 
mance index J and the gradient of the function with respect to 

each variable, i.e., allowable feedback k.,, is provided in 

^D 

the matrix G. A descent technique based on the gradients such 
as FI etcher- Reeves (given in Appendix F) is now used to 
minimise J and to improve the feedback matrix E. This E 
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matrix is then used to start the procedure and the procedure 
is repeated till the condition for optimality is achieved i.e., 
g. . = 0 for each allowahle i and j . 

JL J 

6.4 DESIGU OP GCFSTilAIEED CONTROLLER FOR THE 
PRESSURISED HEAT? WATER REACTOR 

The method described in the preceeding section is now 
applied to the pressurised heavy water reactor to develop an 
optimal controller of specified configuration. The reactor, 
the dynamics of which is represented by Eq.(4.7), has three 
control variables. These are the reactivity u^, the area of 
the feedwater valve Ug and the area of the bleed or feed valve 
for the primary circuit heavy water u^. State variables of 
primary interest for efficient opera.tion of the reactor are 
the following 

i) average temperature of the coolant in the reactor x^ 

ii) water level in the steam drum x^^g 

iii) pressure of the primary coolant circuit x^^g and 

iv) steam drum pressure x^^ 

The design and operation of the reactor requires that the 
first three of the above four variables are to be held at 
their initial values. These variables are easily measurable 
and are readily available for control purposes. Therefore 
the three control variables must be generated from these 
four variables. To achieve this, the following control 
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configuration, wiiich is similar to conventional analog 
controller is postulated. 



t 


^1 

= aj^ + a 2 ^ ds 

( 6 . 14 ) 


t 


^2 

\ ^16 ^2 ^16*^® 

( 6 . 15 ) 

^3 

°1 ^19 

(6.16) 


It is now required to find optimal values for the coefficients 
^ 1 * ^ 2 * ^ 1 ^ ^2 *^ 1 * 

The controller represented by Eqs.(6.14), (6,15) and 
(6.16) is of the proportional - integral type. In Chapter 5, 
the proportional- integral controller has been designed by 
miOdifying the linear regulator theory. Following the method 
outlined in that chapter, the dynamics of the pressurised 
heavy water reactor ca.n be represented by the equations (5.16), 
(5.18) and (5*19) which are repeated below for convenience 

Reactor dynamics 

i = 4 z + B u (5.16) 


[4 

0 , 

\^i 

= i 

; B = 


Ic 

0 


L“ 


“1 


where = (x^, y^) 

z® = (^4’ ^16’ ^19) 

= u ■ 


and 


V 
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Performance Index: J = i ' (5.^ £ z) (5:18) 

and the control law is v = Zz (5.19) 

The proposed controller, represented by Eqs.(6.14),- 
(6.15) and (6,16) is transformed into the system variables* 
and V and the control variables can be expressed as 


*^1 ^ 4 . ^2 ^4 

+b2 

Uj = Vj = 



Therefore, the proposed controller when expressed as v = Kz 

has a specified structure with the elements k, >, > k^ 

JL ^4 

^2 16^ ^2 21^ ^5 19 allowed to take non-zero values. 

The weighting matrices Q and R for the performance 
index in Bq.(5»18) are chosen as 

Q = diag (0, 0, 0, 0,, 0, 0, 0, 0, 0, 0, 0, 0 , 

0, 0, 0, 0, 0, c, 0, i,ai,o.i) (6.17) 

and R = diag (1,1,1) (6.18) 


In order to select the optimal values for the non-zero elements 
of the feed ba.ck matrix, the procedure discussed in Section 6.3 
is now applied to the system represented by Eqs. (5.16) and (5.18) 
with the weighting matrices given by Eqs. (6,17) and (6.18), 

The minimisation of performance inde± J is done with the use 
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of Pletciier-Reeves metiiod (given in Appendix P) . and the 
optimum values obtained for the non-zero elements of the feedback 
matrix F arc given by 


1,4 

^ 2,16 


-1.114 

-1.0535 


k 


U,20 

2,21 


-0.7997 

- 0.00926 


and kj = -19.985 


Since ^1,20* ^2 16* ^2 21* ^3 pg spend to a^^, a^, 

^1* ^2* °1 ^sspectively, the control scheme specified by 
Bqs.(6.14), (6.15) and (6.16) can be written as 




= -1.114 


4 


- 0.7997 


t 

I 

0 


U2 = -1.0533 - 0.00926 


u^ = -19.985 


t 

J 

o 


^4 as 

(6.19) 


(6.20) 


(6.21) 


These equations constitute the optimal controller of the 
specified structure. 


The performance of the reactor with this controller 
providing the feedback control is obtained for a -5 percent 
step change in the area of the steam valve (Fig.6.l). The 
dynamic behaviour of the system with this control law compares 
well with the behaviour of the system with complete state 
feedback . (Fig. 5.4) • However it can be seen that the transient 

responses of the average temperature of the coolant in the 

reactor ATI, are different for the two control laws. Also 
ac 
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tile magnitudes of required control i.e. values of 6k ^ and 

tlian tliose obtained with complete state 

feedback. 

6.5 MCDIPICATIOH TO THE WATER EEVEI COFTROL 

A modification to the level control is envisaged. In 
Eq.(6.15) the feed water- flow rate . is assumed to be proper— 
tiona.1 to the error in the water level and to its integral with 
respect to time. This control scheme may be adequate for 

OWtn, pttvktdS 

normal fluctuations in the water level^ But a sudden decrease 

in drum pressure for example, will result in an initial rise 

in ihe level, owing to an increased discharge rate from the 

riser and boiling legs of the steam generator. This will initiate 

a control action so as to reduce the feed water flow, whereas the 

correct control action should be to increase the flow. To 

circumvent this tendency^^, the error in the steam flow rate 

is added to earlier control, namely Eq. (6.15). Hence 

t 

u^ = x^g + b^ I ^16 "^^'s (6.22) 

Since w^ = 0.1776 x^^ (4# 5 ) 

u^ can be written as ' ‘ , 

'^2 “ h ^16 + ^2 1^16 h .5 * "1 


With this modification and using the values of b^ and b 2 


obtained from Sec. (6.4), u^ is written as . 


t 

, ; 

3 

U2 = -1.0533 - 0.00926 / ' x^gds 

+0.1776 x^^ + w^ 


(6.24) 
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Closed loop responses of the system, for a —5 percent 
step change in the area of steam valve w^, with the control 
action provided hy Bqs. (6,19) » (^. 21 ) and (6,24) are chtoined andare 
shown as dotted lines in Big. 6,1, It is seen that the above 
modification has not reduced the erroneous initial movement 
of the feed water control valve U 2 . This may be due to the 
heavy weightage provided in the matrix £ for the level. It 
is expected that a reduction in the weighting factor for level, 
namely Q.21’21* allow variations in the level and will 

reduce the unnecessary movement ' of the feed water valve. The 
weighting factor ^21 2 I suitably selected such that 

the- initial movement of the control valve is governed by the 
steam flow rate and the level error signal is used to maintain 
the drum level at the correct value over long periods of time, 

6.6 PERBCRMCE CP THE NOHLIHEAR MCEEL CP- THE , PHWR 

In Sec. 6, 5 » a controller of specified structure which 
is quite similar to an analog controller is designed by 
applying the optimal control theory to the linear model of the 
pressurised heavy water reactor. In order to stiady the 
behaviour of the plant with this controller, the controller is 
incorporated in the dynamics of the plant which is represented 
by the nonlinear model (of chapter 5)* The performance of the 
system for a -5 percent step in the area of steam valve is 
obtained and is now discussed. 
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6.6*1 Controlled Response for a —5 percent step in Area of 
Steam Yalve 

The controller obtained in Sec. 6. 5 is represented by 
E qs.(6.19)> (6.21) and (6.24). These equations are repeated 
below for convenience. 


t 


It 

-1.114 - - 

0.799? f ■ 
0 

x^ ds 

(6.19) 

II 

CM 

.irl.0533 x^g- 

0.00926 

I ds 



+ 

0.1776 x^5 

+ w^ 

(6.24) 

5 

-19.985 



(6.21) 


Substitpiting the physical plant variables for the state and 

control variables, these equations can be written as 

t 

= 6k^ = -1.114 At^^ - C.7997 ^ AT^^ da (6.25) 

t 

^2 " ^'^m " -1-0533 Ay - 0.00926 l 'Ay ds 

+ 0.1776 A p +A (6.26) 

u~ = A a^ = -19.985 Ap (6.27) 

3-D 

Since ■ ax represents the deviation of the variable z from 
the initial steady state value, A x .can be written as 

Ax = x(t) - x(o) 

Further Aa^ and ao^ represent the ' deviations as ' perGentages; 
therefore Aa^ can be written as - 

Aa^^ ” [aj^(t) - ajj^(o)] x 100. 0 
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Hence the actual valve opening a (t) becomes 

X 0,01 

Similarly a^^Ct) = ajj(o) + x 0.01 

snd a„(t) = a^(o) + AWt x 0.01 

S S X 

The initial values of theplant variables are given below 

T^^(O) = 52C.68°P 

y(0) = 0.0 ft 

p(0) = 568.9 psia 

P(0) = 1250.3 psia 

a^(0) = 1.0 • 

a^(0) = 0.0 

ag(C) = 1.0 

In the linear model the level y has been expressed in the unit 
— 2 

of 10 inch where a,s in the nonlinear model it has been expressed 

in feet. Making appropriate corrections in y for the change 

of units between the two models and substituting the initial 

values for plant variables, the controller given by Eqs. (6.25), 

(6,26) and (6.27) can be represent:ed as 

t 

6k = -1.114e - 0.7997 Teds (6.28) 

G o 

where e = - 520.68 

t 

a^t)= 12 X (-1.0553 y - 0.00926 I 

+ 0.001776 (p-568.9) + 

‘ D 


(6.29) 
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a^(t) = -0.19985(1-1250.3) ( 6 . 5 c) 

Thus Eqs. (6.28), (6.29) and ( 6 . 3 O) represent the analog 
controller scheme, as applicable to the non-linear model. 

This analog controller is included to the dynamics 

of the nonlinear model which is. represented by Bq.(3.17). 

The response of the system for a -5 percent step in the area 

of steam valve i.e. a^ = 0. 95 is obtained (Pig. 6 . 2). It can 

b® seen from these responses that there is no change in the 

values of neutron power n’ and the average temperature of the 

coolant in the reactor T. . for the first seven seconds. This 

is due to the transport time lag of the coolant between the 

steam generator and the reactor. At the end of 7 seconds, T„„ 

increases sharply due to the increase in the inlet temperature 

of the coolant to the reactor; this introduces a steep decrease 

in the reactivity 6 k: and steep variations in the area of the 

c 

valTSfe opening a^. 

6,6.2 Improvement of the Responses by the Peed back of ' 

Steam Drum Pressure 

The responses obtained in the proceeding section can 
however be improved by modifying the reactivity control. The 
sharp changes in T , and 6 k are mainly due to the 
transport delay between the steam generator and the reactor . 

It is expected that the effect of the transport delay can be 
reduced by including steam drum pressure s». an anticipatory 
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signal in the reactivity control law Eq.(6.28). 

The control law for the linear model incorporating 
steam drum pressure in the reactivity control can be written as 

^ "■ c ^ ^1 ^^ac a 2 ; A T^^ da + a^ A P (6.31) 

^2 “ ~ ^1 + ^2 (6.32) 
U 3 = Aa^ = c^ AP ( 6 . 33 ) 

Optimal values for the constants a^, a 2 , a^, b^i b 2 and c^ 
are obtained for this controller by the method discussed 
in Section 6 . 4 . The cost matrices g and R are the same as 
given by Eqs.(6.17) and (6.18), The optimum values thus 
obtained a3?p 

a^ = -1.383 a2 = -0.2263 a^= -0.17245 

hj|^ — —1.1829 b2 — —0.08803 

and = ■*25 

Therefore the control law for the nonlinear model can be 

written, similar to Eqs.(6.28), (6.29) and ( 6 . 3 O), as 

t 

^^c ~1'583 e -0.2263 Zeds - 0.17245 (p-568.9) ( 6 . 34 ) 

V O. . 


aj^(t)= 12.0 X (-1.1829 y - O.O 88 O 3 ^ y d 

' ' ■ ■ ■ ■ 0 ' 

+ 0.001776 (p-568.9)+ a 

, ,3,, , 

ajj(t)= - 0.257 (P-1250.3) 

The response of the system with this contr 
for a -5 percent step in the area of the steam val' 


^3 


signal in the reactivity control law Eq.(6.28). 


03ie control law for the linear model incorporating 


drum 

pressure 

in 

the 

react; 


= 6k = 


AT 




c 

1 

ac 

2 

U2 

ti 

t> 

II 

h 

Ay 


’=2 

^3 

II 

O 

II 


ap 




t 

I 

o 


t 

; 

o 


(6.31) 

( 6 . 32 ) 

(6.33) 


Optimal values for the constants a^» a 2 > h^i ^2 and C 2 _ 
are obtained for this controller hy the method discussed 
in Section 6.4. The cost matrices ^ and ^ are the same as 
given hy Eqs*(6.17) and (6.18). The optimum values thus 
obtained aaja 


a^ =-1.383 a2= “0.2263 a^ = “0.17245 

b^ = -1.1829 b2 = “0.08803 
and c^ = ^25*7 


Therefore the control law for the nonlinear model can be 
written, similar to Eqs.(6.28), (6*29) and (6.30), as 

■ t 

6k = -1.383 e -0.2263 ^eds - 0.17245 (p-568.9) (6.34) 

^ O t 

a (t)= 12.0 X (-1.1829 y “ 0.08803 ^ y ds) 

m Q 

+ 0.001776 (p-568.9)+ a„ (6.35) 

"O' '■ 

ajj(t)= -0.257 (P-1250.3) (6.36) 

i ' 

The response of the system with this control scheme 
for a -5 percent step in the area of the steam value is 
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shown in Fig. 6. 3* The performance of the system is "better 
than that of the system with the previous control scheme 
and control input a^ is also less. 

6.7 CONCIUSIOES 

In this chapter, optimal analog controllers are 
obtained for the linear model of the press\irised heavy water 
reactor and the dynamic behaviour of the linear model is studied. 
Later these controllers are suitably modified and incorporated 
in the nonlinear model and the performance of the reactor for 
a ~5 percent step change in the area of steam valve is analysed. 



CHAPTER 7 


t summary md suqcestiohs 

7.1 SUMMARY AHDCOHCLUSIOHS 

A ma;ior portion of the present day energy require- 
ments is met by hydroelectric and fossil fuelled thermal 
power stations. However the increasing demand for power 
and the current oil crisis are bound to place more emphasis 
on nuclear power stations. This will require development of 
control system that will increase plant availability and 
prevent outages; that is, there is a need to design control 
schemes that will resiilt in efficient and safe operation of 
the plant. With this aim, controllers are designed for a press- 
urised heavy water reactor. 

The first step to design the' controller is to understand 
the process that is to be controlled. Hence a mathematical 
model of the pressurised heavy water reactor is formulated 
from the conservation of mass, momentum and energy. The 
resulting nonlinear model is simulated on an IBM 7044 digital 
computer to obtain the uncontrolled responses for a variety 
of disturbances. In order to facilitate the design of 
controllers, the nonlinear model is linearised around stead37' 
state operating conditions and the linear model is subsequently 
used for the design of the controller. 



Since there is a tendency to incorporate process 
computers in the operation of power plants, it is feasible to 
use the computers for control purposes also. As a step towards 
this form of control, feedback control schemes can be designed 
based on modem control theory. Hence in Chapter 5, optimal 
control theory is made use of in designing proportional and 
proportional-integral controllers for the linearised model of 
the plant. The behaviour of the system with these controllers 
are discussed and the drawbacks of these controllers are also 
outlined. 

Power stations in general have analog controllers. 
Therefore it would be advantageous to adjust the gains of the 
existing controllers or to modify the control configurations 
so as to improve the performance of the plant, without resorting 
to new control schemes (of chapter 5)* Therefore in Chapter 6, 
optimal values for the gains of an analog proportional -integral 
controller is obtained by using linear regulator theory and 
function minimisation techniques. This controller has been 
designed based on the linear model. Therefore it is modified 
to include the initial operating conditions of the reactof. 

Then it is incorporated in the nonlinear model of the plant. 

The response of the plant with this controller is obtained and 
is found to be satisfactory for the disturbance studied. 
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7.2 SUGGESTIONS POE PUEtHER STUDY 

Some of tDe related problems for further study are: 
i) The fluctuations in the feedwater valve may he eliminated 
or reduced by decreasing the weightage on the level variation 
or by increasing the weightage on valve movements. This will 
allow for fluctuations in the level during the early part of 
the transients with the primary control of level being effected 
by steam flow rate and with the proportional-integral feedback 
of level acting as a readjusting control so as to maintain the 
drum level at the correct value over long periods of time. 


ii) a) In the analysis, the dynamics of the various valves such 
as steam valve, feed water valve and bleed or fsei valve for 
primary coolant may be considered. This will yield a true 
representation of 'the system. 


b) Purther it is realistic to consider the rate of change 
of reactivity, i.e., ^ 6k^ as the control input rather than 
6k^ itself. The analysis then will include the behaviour of 
the control-rod drive mechanism. 


0 c) In this study, feed water flow has been assumed to 
depend only on the area of the valve opening. This can be 
modified as to incorporate the dynamics of the boiler fe-ed 
pump and its fiepoadence on the steam drum pressure. 
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iii) The method of designing optimal analog controller 
can he extended 

a) to adjust the gain of the existing analog controllers and 
h) to achieve a controller that will make the reactor 
follow the turbine load fluctuations, 

iv) Since the primary heat transport system is sensitive to 
energy transfer from the reactor to the steam generator, the 
pressure of this system can he used to provide an anticipatory 
control; is, the primary circuit pressure can he utilised 

to vary the external reactivity. This procedure will reduce 
the dependency of reactivity control on the transport delay 
between the reactor and the steam generator. 

v) lastly, the dynamics of turbine and generator can be 
included in the model to provide a complete understanding of 
the power station. 
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APPENDIX A 


PXIBL DINMIGS 

Major portion of the heat is released within the fuel 
itself. This heat is transferred to the coolant through heat 
conduction in the fuel and the canning (clad). The heat 
transfer from the fuel to the canning is a mixture of convection 
and direct metallic contact. 

The following assimiptions are made in evaluating the 
average temperature of the fuel, 

1. All fuel rods in the core are assumed to have the 
same power. Power generation in the fuel rod is 
assumed to he constant (i,e. independent of radius), 

2. Temperature in the fuel depends only on radial distance 
from the centre line of the fuel rod, 

3. Thermal conductivity of fuel is described by*** 

k = a/(b+T) 

4. The heat transfer from fuel surface to outer canning 
surface is described by one heat transfer coefficient. 

To evaluate++ the temperature distribution in the fuel 
rod, the fuel rod is divided into a number (m) of concentric 
cylindrical shells, each of equal volume. iEhe radius of inner 

+ I. Devoid, ’A Study of the Temperature Distribution in UO 2 
Reactor Fuel Elements’, Report AE-318, Aktiebolaget 
Atomenergi, 1968, . 

-H- Power System Simulator - Description of Models, IBM Manual, 

GE 20-0387'-0, 


1 


A-2 

bo-undary, of shell no. i is then 

'{ — i = 1, 2, , m 
where = radius of fuel rod, 

let tJe "tile average temperature in shell number i. 
Then the temperature gradient on the inner boundary of shell i 
can be written as 


M I 

dr ir=r . 




i * 2,3, ... ,m 


= 0 at i = 1 

The heat balance in shell no. i, then, is 


itr^ dT^ 

pc ■— *- v - x— 

m dt 


+ a-rtk. T A r. 
m 1-1 1 1 


T. .-T. 
1-1 1 




- 27i:kj_ I 


T -T. 

•^i ^itg 


where Q = heat generated per fuel rod 
a = length of fuel rod 
k^ = thermal conductivity of shell umber i 
p = density and 
c = specific heat of fuel 
The above equation can be written as 


i — 2,3,... ,m— 1 


dT. 


dt 


^ ni ^m *** ‘^'f ^^i-1 \-l^'^i-l‘'\^ ^ ^i^i ^^i”®i+l^^ ^ 


i = 1, 2 , . , ,m-l 



wheye 


[^ = 0 % 




= 


fl+i t fl^ 


X ^ lj^2«»t 


^0 * 0 


Fow ■th.e outer most sliell 




dt H. ‘■m ^ f "^m-1 ‘"m-l 2 ""gap^^fs -"co 




U.2) 


■where iigg^p ~ hea't transfer coefficient "between the fuel 

surface and clad outer surface 
T^g - temperature at fuel surface 

- temperat-ure at clad outer surface 
On the boundary we have the following eq.uation 


'■ fs co^ gap m 


^m”^m*fl 


f A,/m+l /m-^> 


U-3) 


Here is "tiis temperature in an auxiliary shell outside ' 

the fuel, 
ilso we have, 

- Tfm^ 

fs m ' m m+1' " " . (i*.4) 

Tm+1 - Y m-1 

Prom Eqs., (A, 3) and (A,4), is solved and introduced in 
Eq, (a, 2), This leads to 


dT„ 


V ^ V-l ’’b-1 ( Vi-\) 


1 . 



A-4 


~ A. 


f ^ m-1 t 

(1-^ ) . - 


gap 


m-l 


(A.5) 


Equations (A.l) and (A. 5) are used to calculate 
i=l,2,,..m with a constant value for T 


ature of the fuel T 


af* 


The average temper- 

can he evaluated as 


1 ^ 

L-f = i- ^ T. 

Six HI ^ 


(A. 6) 


To calculate the fuel temperature distribution the 
number of sections (m) was increased until no significant 
improvement was observed in the calculated values. A value of 
6 was used for m and adjusted so that the average 

temperature matches with the design value at full power. 

Data used: 

h »= 370 W/ft.^ op 

gap 

k = 46x16.9/(480 +(5/9XT-32))^< 2500°F 

= 16.9x0.024 T> 2500 °P 

k in ¥/ft. and T in °F 
r^ = 0.0235 ft. 

p = 660 Ib/ft^ 

c = 59.2 WS/lb 

Q = 7Q00 W/ft. of fuel rod. 

^co = 537 


At full power 


V = 

= 1110 


(calculated based on 6 segments) 
°P (design value) 



Prom the steady state fiiLl power operation, the 
temperature response for a +10^ step in power is obtained. 

The response of the average temperature is shown in Pig.A.l. 
Simplified fuel model: 

The above representation of fuel dynamics requires m 
ordinary differential equations with a low limit for integration 
step about 0,015 sec. A much simpler representation of the fuel 
dynamics is therefore required. Here the average fuel temperature 
Taf, is calculated from the following equation 

dT^f 

"dt ' % "* ^^af -^co^ (A. 7) 

where Aj^ = Heat transfer area at clad outer surface 
= 2%Z-x outer radius of clad. 


Here H^ is an artificial heat transfer coefficient 
which is assumed to be a constant and is calculated with the 
condition that Eqs. (A. 6) and (A. 7) give the same value for 
T^f at full power steady state operating conditions. Hence 
H^ becomes equal to 83 W/ft^ •?, was adjusted such that 

the response of the simplified model compares favourably with 
that of detailed analfsis for a + 105 step in power. The 
responses are shown ih Fig, A.l with '^^=5.0 seconds. 


Pigure A.l shows that the two responses agree well 
during the early part of the transients, later settling to 
different steady state values. This is because, in the detailed 
analysis the thermal conductivity is treated as a function of 
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APPENDIX B 


HEAT EXCHANGER DYNAMICS 

Eor purposes of analysis, the heat exchanger can he 
divided into three portions: the recirculating leg performing 
the role of riser of the natural circulation loop and the 
boiling and preheat legs of the feedwater flow. For the recir- 
culating leg, the following assumptions are made: 


1. Temperature is dependent only on axial distance. 

2. Density variation of the primary coolant (heavy water) 
with temperature is neglected. 

3 . All the primary coolant pumps are in operation at 
the rated speed. The velocity primary coolant in 
the heat exchanger is assumed to be constant i.e., it 
does not vary with time and with radial distance. 

4 . For the secondary coolant (water), nucleate boiling 
is assumed and the heat transfer from riser walls to 


water is taken -fco be proportional to the fourth power 
of the temperature difference between the wall and 
saturation temperature. / 

Heat flux (in Btu/hr ft^) = constant x exp(p/225) 


^^wall^^sat^ 



where 



c - specific heat 

- overall heat transfer coefficient between 
the heavy water and. tube material* based 
tube outer surface area 

- boiling heat transfer coefficient 

(10^60“^) exp (p/225) 

h - enthalpy of secondary, coolant 

w^ - mass flow rate of secondary coolant in the 
riser 

u - velocity of primary coolant 

s outer perimeter of tube 

z - axial distance (along the riser) 

Using Eq. (B,l), Eq. (B.4) can be written as 

^2P fl + ’'r 3 !- = =2 - Vt) 

Equations (B.2), (B»3) and (B.5) predict the +emperature and 
enthalpy distributions in the riser. 

This set of partial differential equations is converted 
into a number of ordinary differential equations by eliminating 
the partial derivatives by finite (backward difference) 
difference technique. The heat exchanger is therefore 
divided into n sections, not necessarily of equal lengths, 
and the energy balance equations are written for the i^^ 
section of length Jl^^, by approximating 
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0?h.en, the i section can be represented by the following 


equations : 


(pc 


i T.-T. 

" + ^i^hw ^ " 

1 


_irl = H s (T .-T.) 
H. 1 ^ W1 




* dt 


i + Wrl = Ha s (^r^sat) 


■(B.5a) 


The subscript i refers to the properties at the i^ 


section. 


The mass flow rate w^^ has to be calculated by solving the 
equation for momentum balance for the entire natural circu 


lation loop. But in this analysis, w^^ has been assumed to be 
a constant. 


The above set of 311 aqiaations are solved for steady state 
distributions at the full load operating conditions (design 
values), by making the time derivatives equal to zero. The 
resultant algebraic equations are solved with the following 
values; 

Inlet temperature of heavy water = T^ = 560 °F 

pressure p = 583 psia 

w^ = 41*5 lb /sec 

The value of n is increased ■until no significant improvement 
is observed in the calculated distribution, A value of 5 for 
n and = 2, 2, 4, 6, 6 ft. for 1 = 1,2, ...5 are used. 
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From these steady state values, the single energy storage 
element approximation for the riser leg represented hy outlet 
properties* is made. The heat transfer coefficient and 
the constant in are calculated such that both the n-^segment 
approximation and the single energy-storage element approximation 
yield the same outlet properties for the primary and secondary 
coolants. These values of and are used to represent the 
riser leg and its dynamic behaviour is hence represented by 
the following equations: 


0.43 - 1*61x1.30 (T,-T .) 


11 2 
JS.61 


2 *w2' 

5 1.61 


^w2 = ^ Ott ^ Of 


h = exp(p/2S5) (\2-l3atl+ "r < 

/(13.1 P^) (B.6) 

where Tg* ’^■w2* ^r temperatures of heavy water and 

of wall of the tube and the enthalpy of the secondary coolant ar 
the outlet of the riser and T^ refers to the temperature of heavy 
water at inlet to the riser leg. 


-3 


Transient responses of the single lumped system and the 
n-segment system are compared for various disturbances. I'he 
disturbances considered are: 

1. a + 10 °F step in the temperature of heavy water at d 
irO-et to the riser, with the drum pressure at 583 psia,and 

2. a +17 psia step in the drum pressure with the inlet 

temperature of heavy water at 560 The responses 
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are shoim in Figs. B.,1 and B.2 and the agreement is 
quite good., Hence the dynamic behaviour of the riser 
leg will he represented hyEq, (B. 6 ). 

A similar study is made for the boiling leg. (Ihe 

operating conditions are taken as: pressure = 583 psia, 

temperature of the heavy water at inlet = 510 enthalpy of 

the feed water at inlet == h„, and mass flow rate of the feed 

f 

water *= 8.7 Ib/sec) . Ihe responses of single ,l\imped system and 
5 -segment approximations are in good agreement for disturbances 
on pressure and inlet temperatfure of the heavy water and hence 
the boiling leg can be represented by: 

= 0.577 ( 12 -^ 3 ) - 1.32x1.32 

m — ii.22xI^7A ( m ) _ Sxl^^x ext) (n/225)(T 

0.85 ^ 3 w 3 '^. 02 J.U X exp 

h (p/225) 

“ "'m ^ ^ 0.102/ (B.7) 

r 

where I 3 , 1 .^ 3 , h^ refer to the ' temperature of heavy water at 
outlet, the temperature of the wall at the outlet end of the 
boiling leg and enthalpy of feed water at it g outlet end 
respectively. ^2 tcef ® 3 rs to temperature of heavy water entering 
the boiling leg and is equal to the temperature of heavy water 
leecving the riser leg. Similarly h^ refers to the enthalpy 
of the feed water entering the boiling leg and w^ refers t© 
the mass flow rate of the feed water. 




iu 

i^■■ 


Time , sec. 

Response Ror a 4, fO'^F Step In The Iniet Temperoture Of 
Heavy Water 


Time, sec 
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She preheat leg of the heat exchanger is analysed in 
a similar manner. She number of sections considered, is ten 
and tile energy balance equa'^ion for "^he feed water is modified 
to allow boiling or convective heat transfer depending on the 
enthalpy of the feed water. Dynamic behaviour of the preheat 
leg approximated, as a single lumped section, can be represented 
by the following equations: 


i. / 


s 


w4 


= 1.76 (lys^) - 0.85x0.377 
= 1.74x0.377 - 5.0x1.19 

5t+ 


T H-+S . 
sat ,mi 


) 


■1' .AX' 


[■ 171x1 .19 (X^ 4 - 


2 ) - 

/(1. 39 pf) ■ '(B.8) 


Here 5/ refers to the temperature of the heavy water at 

•i. 

the exit of the preheat leg, refers to a hypothe'^ical 

temperature of tube wall, h refers to the enthalpy of feed 

ir 

water at exit of preheat leg. refers to temperature of 

heavy water at inlet to the preheat leg and refers to 
the temperature of feed water at inlet to the preheat leg. 
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STExtfl DRUM , 

Steam drum is idealised as a cylindrical vessel with, 
flow rates as shown in Eig. C.l 

Symbols: 

h = enthalpy 

V = specific volume 
M = mass in the drum 
w = mass flow rate 
p .= drum pressure 

Tgat ~ saturation temperature at pressure p 
p = density 

Subscripts: 

f = saturated liquid 

g = saturated vapour 

mo = feed water from boiling leg 

r = riser 

d = downcomer 

. (dot) indicates differentiation with respect to time t. 

To evaluate the dynamic behaviour of the steam drum, the 
following assumptions are made; 

1. 1. There is no spatial variation of temperature in the 

drum and the temperatures of liquid and vapour phase 
are assumed to be eq ual to saturation temperature. 
This assumption is valid since the feed water and the 




^lii: 








t7m\s 
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riser flow are at saturation temperature and if tlie 
steam separation is perfect. 

2. liquid level changes due to huhble formation in the 
liquid a.re neglected. Further, variation of the liquid 
surface area with changes in the level is neglected,,., 

i.e., the drum walls are considered to be plane and 
parallel for small changes in level and 

3. Heat capacities of the drum internals and the drum 
material are neglected. 

The first assumption implies thermodynamic equilibrium 
between the phases during transients. It is possible to relax 
this assumption and formulate a model which includes mass 
transport at the liquid surface"^. But in the absence of 
experimental results on steam drum alone to verify those models > 
it is considered acceptable to assume thermodynamic equilibrium. 


Mass and energy balances for the contents of the drum 
yield the following three equations for the three variables 


^ . 

dt ’ 


dt 

Mass balance: 


dM. dM 

and 


§M 

dt 


M = M^ + M_ = w^^ -• w„ + w„ 
X g mo s r 


- w- 


;c.i) 


1. K.l. Chien et al., 'Dynamic Analysis of a Boiler', Trans. 

ASME, Vol.80, pp. 1809-1819, 1958 

2. P.Z, M'Bherson et al., 'Dynamic Analysis of a Nuclear Boiler', 
Proc. IMS, Vol. 180, pp. 417-427, 1965 

3. P.E. M'Pherson, 'A Mathematical Model for Steam Drums', 

UKAEA Report ABEW-R 366, 1966 
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Energy balance: 
d 


dt (“A + W'> - It (5^) = ’"mo\ + ''t\ ~ "d^ ■• '‘b\ 

neglecting the second term on the left side, the above 


equation ca.n be written as 

dh^ 


Vf + + (“f dF ^ “g-df") H 




d f 


s g 


( 0 . 2 ) 


Volume of tho driim is constant and this is expressed as 
V = + ^ 0 ^ 0 - - constant 

Hence 

H = Vf + Vg + (“f 


dv^j dv^ , 

-f-. . M — ^ 

dp ^g dp dt 


0 


(C.3) 


Equations (C.l), (C.2) and (C.3) are solved, after 

dVf dh„ 

neglecting terms involving -x — and M v— “ for the three 

^ up g up 

unknowns viz., p, M„ 

I g 

= B ' /C ' 


and 


dt 

dM^ 

dt 


D>/E’ 


(C.4) 

(0-5) 


where B' = w^ h^ - w^ (h^ - h^) 

"■ ’'no <’'f« ■ ’'f *’fg^ ■ '’'r - ’'d) 

0’ = Mg {hjg - P) 

= ’'no^VN’ ■ ’'d 

- [ (’'no-’'s) + K-’'d> 3 V' 
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E' 

E’ 


and a. 


= 

= Mf aj^/M aj 


8*1 ““ 


dp 

d'^n. 

s 

dp 


Further , ' 

= p V 
f f 

where = volume of liquid in drum. 

Measuring the liquid level y from the centre of the drum, 
can he written for small variations of y about the centre 
of the drum, as 

Vf = ( “. + y.d) £ 


where£ = longitudinal length of the drum 

2 

_J_ TT/^ ^ £ , 

f 


Hence 


M. 


(S|l + yd) I P. 


and M 
dM 


8 

= (5|1 - yd) £ 


f 

dt 


p 

f dt 


This yields dy/dt as 


dt 


p . £d 


g 


E' 


■( 0 . 6 ) 


and D* and E' are given in Eq. (C.5). 


The above Eqs. (C.4) and (C.6) involve the flow rates 
from the riser and boiling legs and the flow rate to the 
downcomer. Following section develops expressions for these 
flow rates. 
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lJat*ural circulation loop: 

Steam drum, downcomer and riser fwiiicH’ Is the .shell of 
one leg of U-tuhe type heat exchanger) form the natural 


I 


dil^dulatibn ibopi Mhis flo^ rale in Ihis can Ife etulied ly 


considering the momentum "balance for this circuit'*'. 
The following assumptions are made: 


1, Ho "boiling takes place in the downcomer and the fluid 
in the downcomer pipes is saturated liquid. 

2» Vapour and liquid velocities in the riser are equal. 
Also thermodynamic equili"brium is assumed to exist 
between the phases. 

5. Kinetic , energy of the riser fluid entering the drum 
is dissipated in turbulence and the liquid in the drum 
is stationary (has no velocity towards the downcomer 
pipes). 

4. Data for pressure drop due to internals of drum (like 
steam separator) are not available and hence that is 
not considered. 

Momentum balance for downcomer pipe; 

Available pressure loss due + gravitational + entrance 
pressure = to friction head loss 

head 

+ inertia loss + head available due to kinetic 

energy at exit of downcomer. 


K.L. Chien et al . , ’Dynamic Analysis of a Boiler ' , 
Trans. ASMS, V0I.8O, pp. 1809-18l9> 1958 
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(i.e) 


(P“Pi 


= f , (- 
d 


d ^d 


;) ^d ~ ^ f ^ d 


d f 


^d g a 


d 

2 

d P f 


( 0 . 1 ) 


where = frictional coefficient for downcomer 

= entrance loss coefficient for downcomer 
= length of downcomer pipe 
a^ = area (flow area) of downcomer 

d^ = diameter of downcomer 

p = drum pressure 

p-j^ = pressure (static) at the hottom of downcomer pipe. 
Momentum "balance for riser leg can "be written as 


p p = f. (-|Ji£~. 

^ a^ d^ g 


\ k 

--)+ ' "idz +-£--5^^- 

- o ' g p. 


Ar . ^ 

w + ( ? 

a^ g r ^ p g 


"d 'i 

Pf ^ 


(C.8) 


where k_ = exit loss coefficient for riser 


f^ = frictional loss coefficient for two phase flow 
in the riser 

= 1+2400 (x/p)*^'^^ p in kg/cm^"^, x - quality of flow 

= frictional loss coefficient for single phase flow 

Figure C.2 shows the variation of the density of riser 

it 

fluid with axial distance. In this analysis ^ has been 


+ K.M. Becker et al., Aktiebolaget Atomenergi Report AE-69, 
70 and 85, 1962 
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approximated to where p ^ is density at outlet of the 

riser. Ihe effect of this approximation is to increase the 
available gravitational head resulting in an overestimation 
of the flow rate. In Apprndix B, it ha,s been assumed (assum- 
ption no, 4) that tjie heat flux to the riser fluid is dependent 
only on pressure and wall temperature and is independent of 
the flow rate of'the fwo-phase mixture, Therefore the enhanced 
flow rate does not affect energy transfer, but may influence 
transients in the liquid level of the drum. 

Mass balance for the fluid in the riser treated as 



and Dg = Vjg/hfg 

(3. w 

Prom Eq. (C.IO) , . can be written as, neglecting the 
second derivatives, 
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^ - z (c.ll) 

where 

^ th P + >=2 - 2Vr i\-^2i) V 

^ = -I? = - 0.1681X10^ 

®2 “ ■ = °-22-T8 

Combining Eqs . (C.7)> (C.8), (C.IO) and (C.ll) , to eliminate 
and p and substituting values for physical dimensions and 
constants, the flow rate ,w^ can be written as 

w^ 

Wd = [20 (P^-P^) - w^ - ^(1.036 + 0.274 X 

X . 96 P p 

(l+2400(^) )p— }+0.95z] /5-51 

( 0 . 12 ) 

where z is defined in Eq, (C.ll) 
and w^ is given by Eq. (C.IO). 

The flow rate entering the steam drum from the boiling 
leg is obtained by mass balance for the boiling leg and is, 

similar to Bq. (C.IO) written as 

2 , * 

''mo = ''mi + \ ft <'’3 5 + ’'4 

where 

^ dh~ 

”3 = 35 ("fg/J'fg) ■ v? 

and = flow rate entering the boiling leg which is assumed 
by neglecting the density changes in the preheat leg, to be 
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equal to the mass flow entering the preheat leg. 

Equations (C. 4), (C.ll) and (C.12) are combined to 
yield, after rearranging the terms 


dp ^ 
dt 

(A'-B'-C’-D'b^ ^)/(E'+E'b^+G'b^) 

(0.14) 

where A* = 

’'s hg 


B' = 

(■^ig h - '^fg h + *^g) («ml + 

\ Pb ^4 ^b) 

C = 

(h^-h^) ''fglvV r h h) 


D' = 

v^ h _ V p ^ 
f fg r^ r 


E' = 

M (dv /dp) (h„ -v„ E'O 
g g^ ' fg fg 


E>' = 

dh_ dv 

dp ''f / dp ^ 


p» = 

(^fg h - ^fg h ’"fg) \ 


&• = 

((V*"f) "'fg + "'f ’"fg) ’’r 


and b^,b2,b^ 

and b^ are defined in Eqs . (0*10) 

and (0.13) 


Mass and energy balance for the riser and boiling legs, 
represented as a single lumped energy storage systems, yield 
(Appendix B) 

^ ^ (p/225) 

/(13.1 ^ (0.15) 

and 

=[ fcif == 0-00®25 exp (p/225) (\r"sat)^-’'ai 

(h^-hp)] X 0.102/ '(0.16) 

Similarly, the drum liquid level can be represented by 

= (A’+B’-C‘-D' )/(E’xE' ) 



A-27 


where 


A' 


’'mo 

w ^ 
mo 

=: 

’'mi + "b ("4 + '>3 

B’ 

=r 


w 

r 


’'d + fr ('’a ii- + I’l P) 

0’ 


"d ^fg 

D’ 

= 


E’ 

= 




dh„ dv 

E” 


dT d?^> 

F’ 

= 

Ad 


and and are defined in Eqs . (C.IO) and (0.15). 

Equations (C. 12), (C.15), (C.16), (C. 14) and (C.IV) 
describe the dynamic behaviour of the steam drum and natural 
circulation loop. 

Further, thermodynamic properties needed for the analysis 
are assumed as a linear function of pressure over the range 
of 450-?700 psia and are given by 
h^ = 540.24 + 0.2178 p 

h = 12fl - 0.01556 p 

o 

^sat= 372.847 + 0.1878 p 

v„ = 1.6782 - 0.0014953 P 

Vj = 0.01785 + 0.5829 x lo'^p 

Figure 0.5 shows a comparison of the values calculated using 

^ * 4 ”' 

the above equations "with those provided by Keenan and Keyes. 

+ J.H. Keenan and F.G. Keyes, ’Thermodynamic Properties of 
Steam’, John Viley and Sons, Inc., 1936. 
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IPPEHDIZ D 


COmOLLABILITY 

A system is said to "be completely state controllable 
if it is possible to transfer any given initial state 
x(t^) to any desired final state x(t^) in a finite time 
interval t^ < t < t^ with an unconstrained control vector u(t). 

Por a linear time invariant system given by 
X = A X + Bu 

X is a (nxl)' vector of state variables 
u is a ■ (mxl)" vector of control variables 
A, B are matrices of appropriate size 
the complete state controllability implies"^ that the composit 
¥-matrix 

¥ = (BtABlA^B:** I is of rank n. 


+ K. Ogata, ’State Space Analysis of Control System' 
Prentice Hall, 1967 



APPENDIX, E 

LINEAR EBG-ULAIQR IHEORY 

In the linear regulator problem'*', vre have a linear 
differential system 

X == ^(t)x + B(t)u , x(tQ) = Sq (E.l) 

where x is a vector of n variables, u is a vector of 
m variables, A and B are matrices of appropriate sizes.- 
and wish to find the control which minimizes the cost function 
(for t- fixed) 

tf. 

J » I x^(t^) S x(t^) + I ^ ^x^(t) £ x(t) + 

^o 

(t) R 2i(t)) dt (E.2) 

The matrices Q,, R and S are assumed to be symmetric. ¥e 
may obtain the solution to this problem via the Pontryagin's 
maximum principle or the Hamilton-^Jacobi equation. TJsing the 
former mothod, the Hamiltonian is written as 

H [x(t), u(t), X(t), t ] = I xT S + I u''^ £ u + 

T ' T 

X 4 X + X B u (B.5) 

Application of the maximum principle requires that, for an 
optimum control, 

1*1 = 0 = I(t) u(t) + B^{t) X(t) (E.4) 

and = ^*X = _Q(t) x(t) + 4 ^;(t) X(t) (E.5) 

3.S - ■ - 

with the terminal condition 

= S x(t^) (E.6) 

+^A,P. Sage, 'Optimum' Systems Control’, Prentice-Hall, 1968 
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Thus u(t) = -R-^ (t) B^(t) ^ (t) (E,7) 

To convert this to a closed loop control, it is assumed that 
the solution for ^ is similar to Eq,(l 46 ).. ' 

X(t) = P(t) x(t) (E. 8 ) 

Substituting this relation into Eq. (E.l) and (E.7)» 

Eq*(E(l) can be written as 

i = A(t) x(t) - B(t) R“^(t) B^(t) P(t) x(t) (B.9) 

Also from lq» (B. 8 ) and (E.5)»^can be written as 
1 = £ x(t) + P X 

= - £(t) x(t) ~l^(t) P(t) x(t) (E,10) 

By combining ■ Bq, (E.9) and (E.IO), we have 

Ci+B(t) A(t) + A^(t) P(t) - P(t) B(t) R"^(t) B^(t) P(t) 

+ x(t) = 0 (E.ll) 


Since this must hold for all nonzero x(t) , the term multiplying 
x(t) must be zero. Thus the P matrix, which is an n x n 
symmetric matrix and which has n(n+l )/2 different terms, must 
satisfy the matrix Eiccati equation- 

P = -P(t) A(t) - /(t) P(t) + 

P(t) £(t) R“^ (t) B^(t) P{t) ~S(t) (E.12) 

with a terminal condition given by 

P(t^) = - ^ (E.15) 

Thus we may solve the matrix Riccati equation backward in 
time from t^ store the matrix 

K(t) = - R~^(t) B'^(t) p(t) 


(E.14) 
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A closed loop control is then obtained from 
li(t) = E(t) s(t) 

which yields the desired state feedback control law. 


(E.15) 



APPEHDIX P 


PUECTIQN MIEIMISATIOE 

KLetcher-Reeves metliod'*' is one of the methods used for 
minimising a function, where a large amount of local infor- 
mation about the function in terms of the gradient of the 
function being minimised, is utilised. 

The problem of function minimisatiop can be stated as: 

To find the vector x such that the function P(x) is a 
minimum where there are no restrictions on the choice of z. 


The gradient of the function F lies in the direction of 
greatest rate of change of fionction and has that rate of 
change as its magnitude. The gradient, which is a vector and 
denoted by ? F, is 'defined as 


VF = G ~ 


3 P 


3X2 * * * 3 X. 


) 


(P.l) 


n 


With this definition, Fletcher-Reeves method can be stated as; 

Selecting an initial choice of x as x^ , the gradient G^ at 

x^ is calculated. The direction in which the function will 

decrease, is then 

S = - G 
-0-0 

The function F is minimised in this direction such that at 
■ x^?- 3 ^ + a' S (P.2) 

F(x 2 _) <P(Xq + S^) for any a. At the new point the 

gradient G, is calculated and the new direction for function 

•“i #' . ' ' 

minimisation is written as 


+ R.L. Fox, ’Optimisation Methods for Engineering Design’, 
Addison-Wesley Pub. Co., 1971 
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! S I 

and the function minimisation in the direction ^ is then 
carried ont to find a' and the new x. Thus at the end of . 
i^^ iteration 





(I'.4) 


and the procedure is repeated till the required accuiracy is 
obtained. 


a' in Eq. (E.2) is obtained by using the method of 
quadratic interpolation. Considering any vector S and 
the initial point with function value E^, a new vector 
X is defined as 


X = Xq + a S. 

The parameter a is then adjusted such that 
Fi = T'(So) 

^2 — ^^—0 + a S ) < 

and P^ = + 2a S)> P^ 

Approximating the function P(a) by H(a) which is a quadratic 
in a, i.e., H (a) = a + ba + ca , leads to, the minimum of 
H(a) at a' 

4P. - 3Pp - P^ 

where a' = — ^ a 
4^2 - 2q - 2P3 

x' cal Ciliated as x^ + a' is then taken as the minimum 
of F(x)i in the direction S. // 



